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ABSTRACT 

Context. The onset of massive star formation is not well understood because of observational and theoretical difficulties. To find the 
dense and cold clumps where massive star formation can take place, we compiled a sample of high infrared extinction clouds. We 
observed the clumps in these high extinction clouds in the 1.2 mm continuum emission and ammonia with the goals of deriving the 
masses, densities, temperatures, and kinematic distances. 

Aims. We try to understand the star-formation stages of the high extinction clumps by studying their infall and outflow properties, the 
presence of a young stellar object (YSO), and the level of the CO depletion. Are the physical parameters, density, mass, temperature, 
and column density correlated with the star-forming properties? Does the cloud morphology, quantified through the column density 
contrast between the clump and the clouds, have an impact on the evolution of star formation occurring inside it? 
Methods. Star-formation properties, such as infall, outflow, CO depletion, and the presence of YSOs, were derived from a molecular 
line survey performed with the IRAM 30m and the APEX 12m telescopes. 

Results. We find that the HCO*(1-0) transition is the most sensitive for detecting infalling motions. SiO, an outflow tracer, was mostly 
detected toward sources with infall, indicating that infall is accompanied by coUimated outflows. We calculated infall velocities from 
the line profiles and found them to be of the order of 0.3-7 km s"'. The presence of YSOs within a clump depends mostly on the 
clump column density; no indication of YSOs were found below 4 x 10^^ cm"^. 

Conclusions. Star formation is on the verge of beginning in clouds that have a low column density contrast; the infall is not yet 
present in the majority of the clumps. The first signs of ongoing star formation are broadly observed in clouds where the column 
density contrast between the clump and the cloud is higher than two; most clumps show infall and outflow. Finally, we find the most 
evolved clumps in clouds that have a column density contrast higher than three; in many clumps, the infall has already halted, and 
toward most clumps we found indications of YSOs. Hence, the cloud morphology, based on the column density contrast between the 
cloud and the clumps, seems to have a direct connection with the evolutionary stage of the objects forming inside. 
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1. Introduction 

Massive star formation is thought to occur in clumps that are 
deeply embedded in the dense gas and dust of the natal molec- 
ular cloud. In these cold and dense environments, observations 
are often hindered by absorption or high optical depth of the 
emission and the confusion of many objects due to a limited an- 
gular resolution. Molecular transitions in the submm wavelength 
range are particularly useful probes to study these regions, since 
their critical densities are similar to the densities one expects for 
regions in which stars are forming. Even without resolving the 
star-forming clumps, one can find signposts of star formation, 
such as outflows and infall, from the kinematics derived from 
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the shapes and radial velocities of molecular lines. In this paper, 
we try to use the molecular emission observed toward the sample 
of high extinction clouds, whose con t inuum and ammonia obser- 
vations were discussed in iRvgl et all (1201 Oh . hereafter Paper I, to 
study the infall properties of the clumps and other indications for 
the onset of star formation. 

Crucial for the understanding of massive star formation is 
observational evidence for the infall of matter. In low-mass star- 
forming regions (SFRs) such as the Bok globul e B335, infall ha s 
been inferred from the observed line profiles (iZhou et al.lll993h . 
There are but few infall studies for candidate high- mass stars 
dWu et al.l2005l:lFuUer et al.i2005HBeltran et al.l2006h : however, 
their number is curre n tly inc reasing rapidly. A recent study by 
iLopez-Sepulcre et al.l (|2010b. using a sample of infrared dark 
and infrared bright clumps from the literature, found a good 
correlation between the clump dust mass and the outflow mass, 
irrespective of the clump being infrared dark or bright. These 
authors claim that the outflow properties during star formation 
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remain unchanged within a large range of masses and infrared 
properties of the clumps. 

For our infall study of dense clumps, we used not only 
the shapes of the HCO^(l-O) line but also the shapes of the 
HCO+(4-3), CO(3-2), and CO(l-O) emission Hnes. The dif- 
ferent molecular transitions probe different densities and there- 
fore different parts of the clump. Lower-/ CO lines trace the 
kinematics of the lower density material of the clumps. At high 
densities, n > 10'*cm"^, and cold tem peratures, ~15K, th e CO 
molecule depletes from the gas phase dKramer et alj|1999l) . The 
level of CO depletion therefore indicates if the clump is cold 
and dense, and hence young, or if it consists of akeady more 
processed material. Apart from the presence of broad wings in 
lines from CO and other molecules, outflows can additionally 
be traced by an enhanced abundance of SiO, which is a good 
shock tracer ( Schilke et al..,1997,) . We also observed H2CO and 
CH3CN, which trace the more evolved stage of the young stellar 
object (YSO). 

Star formation is also evident by the presence of heating 
sources: a YSO's radiation increases the temperature of the 
dust and gas surrounding it. This is not only visible in infrared 
emission of the dust but also in several molecules like H2CO 
and CH3CN. These molecules are often observed toward hot 
molecular cores, because the elevated temperatures (> 100 K) 
in these regions evaporate the icy grain mantles, increasing the 
gas phase abundances of such species by orders of magnitude 
(Mangum & Wootten 1993; Olmi et al. 1996; van der Tak et al. 
I2OOO ). Additionally, for CH3CN one can estimate rotational tem- 
peratures of the hot phase from the emission in lines arising 
from different K levels that cover a wide range in temperature, 
but are close together in frequency. Using these tracers, we de- 
fine starless clumps as clumps without YSOs and signs of in- 
fall/outflows and prestellar clumps as clumps without YSOs but 
with infall/outflows. 

The higher density components of the clumps were stud- 
ied with N2H+ and H'^CO+; N2H+ is known to be a reli- 
able pr obe of cold gas w ith lower depletion than most other 
species (iTafalla et al.ll2002l) . The hyperfine structure (hfs) of this 
molecule also allowed to estimate its excitation temperature and 
column density. The optically thin H'^CO^(l-O) line was used to 
determine the local standard of rest velocity (vlsr) of the clump. 

We tried to correlate star-formation behavior in the clump 
not only with surface density but also with the morphology 
and physical conditions of the cloud. In Paper I, we selected a 
sample of high extinction clouds based on our extinction maps 
from the first Galactic quadrant and performed 1.2 mm dust 
continuum observations of them. We defined three categories 
of clouds based on their 1.2 mm dust continuum contrast be- 
tween the clump column density and the cloud column density. 
Clouds with a low clump to cloud column density contrast (CAf,,^ ) 
Cn^^ <2 were defined as diffuse clouds and thought to be the 
youngest clouds: they are the coldest and show few signs of star 
formation, such as masers and infrared emission. Clouds with 
higher contrast 2<Ca(„^ <3 were defined as peaked clouds having 
at least one clump (or peak, hence 'peaked' clouds) with a high 
contrast to the cloud emission; these clouds were thought to be 
the following stage. They show slightly higher temperatures and 
wider line widths (larger turbulence). The most evolved clouds 
were the multiply peaked clouds, defined by a still higher con- 
trast Cn^i^ >3 and generally containing several clumps (hence 
'multiply peaked' clouds). Most of the maser emission and in- 
frared sources were found toward these clouds. Even though sev- 
eral clouds, mostly the peaked and multiply peaked ones, are 



elongated or filamentary, this phenomenological feature was not 
taken into account in defining the three cloud categories. 

The high extinction cloud sample that was studied in Paper 
I contained 50 clumps, of which 12 were in diffuse clouds, 19 
in peaked clouds, and 19 in multiply peaked clouds. Almost all 
clumps showed NH3(1,1) emission, indicating the presence of 
high-density gas (Paper I). In the current paper, we report the 
results of a molecular survey on the clumps in high extinction 
clouds. With these results, we try to pinpoint their evolutionary 
stage and test our proposed evolutionary sequence of the clouds. 
We observed all the positions that contained NH3(1,1) emission 
with the IRAM 30m telescope in the H1^CO+(1-0), SiO(2-l), 
HCO+( 1-0), CH3CN(5-4) K^O, 1 , 2, 3, and 4 levels, N2H+( 1-0), 
J8o(2-l), and CO(2-l) molecular transitions. Next, a subsam- 
ple of these sources was observed with the Atacama Pathfinder 
Experiment submillimeter telescope (APEX) in higher J tran- 
sitions: N2H+(3-2), H2CO(4o3-3()4), HCO+(4-3), and CO(3-2). 
The APEX targets included mostly active and evolved sources, 
which are expected to emit strongly in these higher excita- 
tion transitions. In addition, a few diffuse and peaked clouds 
were added to allow comparison between the three categories 
of clouds. The line d ata were interpre ted using the RADEX ra- 
diative transfer code d van der Tak et aU2007.) to arrive at models 
of the physical parameters of the clumps. 

The results of the line observations are compared with 
results from similar surveys toward other samples, e.g. 
the study of line emission toward Infrared Dark Clouds 
(IRDCs) and other massive star-forming regions bv'Rag an et al. 

(2006) , Motteetal. (2007), Pillai et al. (2007), Pirog ov et al. 

(2007) , Purcefletal. (2009), Lot^ez-Sepulcre e t alJ (120 lOl) . and 
Vasvunina et al. (2011); the survey of methanol maser sources 
by Purcell et al. (2006); the h igh-mass protoste liar object s sur- 
vey by F uUer et al.. (.2005h and lBeuflier & S ridharan' flOOf); and 
similar work toward hot molecular core s by Arava et al. (2005]) 
and UCHii regions bv Churchwell et al] ( |T9 92). This comparison 
will allow assessment of the differences in evolutionary stages 
covered by these studies. 

The observations and data reduction are described in Sect. 2. 
Sect. 3 gives the results of the infall study, the derived tempera- 
ture and density estimates, the search for YSO indications, and 
the CO depletion study. These results are interpreted in light of 
an evolutionary sequence of the three classes of clouds in the 
discussion (Sect. 4) and summarized in Sect. 5. 

2. Sample of high extinction clumps and 
observations 

2. 1 . The sample 

The high extinction clouds were selected from c olor-excess 
maps in the 3.6yum^.5yum Spitzer IRAC bands ( iFazio et al] 
2004). The maps cover Galactic longitudes 10° < / < 60° in 
the first quadrant, 295° < / < 350° in the fourth quadrant, 
and -1° < b < l°in Galatic latitude. The method to con- 
struct the extinction maps is described in Paper I. We selected 
these clouds based on a 3.6 yum-4.5 jjm color excess (CE) above 
0.25 mag, which corresponds to a visual extinction 20.45 mag 
by Ay - 81.8 X C£'(3.6p m - 4.5 fim) using the reddening law of 
llndebetouw et al.l (l2005t) (for more details, see Paper I). To guar- 
antee visibility with the Effelsberg 100m and the IRAM 30m 
telescopes, we focused on the clouds in the first Galactic quad- 
rant. In Paper I we analyzed the dust continuum emission and 
the ammonia inversion transitions of the clumps in high extinc- 
tion clouds to derive their physical properties: the clumps have 
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masses between 10 up 400 M©, they are cold, ~16K, and are 
found at distances between 1 and 7kpc, with the majority being 
at 3 kpc. A summary of the properties of the diffuse, peaked, and 
multiply peaked high extinction clouds, based on Paper 1, is pre- 
sented in Table [1] Additionally, we calculated the surface den- 
sity, following Lo pez-Sepulcre et al.l fcOlO) Eq. 1: S = 4M/jtD^, 
where M is the clump mass in grams and D the diameter of the 
clump in cm. 

2.2. IRAM 30m observations 

The spectral line survey of the clumps in high extinction clouds 
was performed with the IRAM 30m telescope in 2007, June. 
We exploited the possibility of the ABCD receivers to observe 
simultaneously at two frequencies, one at ~ 100 GHz and the 
other at ~230 GHz. With two different receiver setups, we ob- 
served a total of seven molecular lines, listed in Table |2l The 
SiO(2-l) and H^^CO^(l-O) transitions were observed simulta- 
neously in the same backend, since they are only separated by 
100 MHz. The higher frequency CO(2-l) and 18o(2-l) transi- 
tions were observed in both setups, which increased the signal- 
to-noise ratio by V2. For the ~ 100 GHz lines, we used several 
VESPA backend settings: one with a (relatively) narrow band- 
width of 40 MHz and a channel spacing of 0.02 MHz, corre- 
sponding to velocity resolution of 0.08km s and two wider 
bandwidths of 120 and 160 MHz using a channel spacing of 
0.04 MHz, corresponding to a velocity resolution of 0. 16 km s ' . 
The wider bandwidths were necessary to allow the SiO(2-l) 
and H'^CO+(1-0) lines to be observed together in the 160 MHz 
bandwidth and to enable the broad, multiple-ZT CH3CN(5^) 
line profile to fit within the 120 MHz band. For the -230 GHz 
CO(2-l) and J8o(2-l) lines, which both have very wide line 
profiles, we used the 1 MHz backend, which offers a bandwidth 
of 5 12 MHz and a spectral resolution of 1 .5 km s"' . Table[2]gives 
an overview of the bandwidth and spectral resolution for each 
transition. We observed in position-switching mode, where the 
off positiorQ was located 800" away. 

During the observations, we performed a pointing check ev- 
ery 1 .5 hours on a nearby quasar or on the Hii regions GlO.6-0.4 
or G34.3+0.2. The pointing was found to be accurate within 4". 
The focus check was usually performed on Jupiter or 3C 273 at 
the beginning of each observing run. The opacity at 230 GHz 
was variable from excellent winter weather conditions (t = 0.1) 
to average summer conditions of t = 0.5. The system tempera- 
tures at 230 GHz ranged from 270-920 K. At 100 GHz the opac- 
ity ranged from 0.04-0.1, and the system temperatures were be- 
tween 87 and 173 K. 

The observed output counts were caUbrated to antenna 
temperatures, T^, usiri g the standard chopper-wheel technique 
dKutner & Ulichl1l98lh The T* temperature is the brightness 
temperature of an equivalent source that fills the entire In ra- 
dians of the forward beam pattern of the telescope. This can be 
converted to a main beam brightness temperature, Tmb, by mul- 
tiplying by the ratio of the forward efficiency, Fss, and the main 
beam efficiency, Beff : 

Tmb = ^ X . (1) 

Oeff 

Table |2]lists the efficiencies for each transition. 



' The off position is an observation of 'blank' sky, which is observed 
to remove the instrumental bandpass structure. 



2.3. APEX observations 

The observations were carried out with the APEX 12m sub- 
millimeter telescope located on the Chajnantor plateau in the 
Atacama desert (Chile) during several runs between 2007 June 
10 and November 2. We used the double sideband receiver 
APEX-2A equipped w ith two fast Fourier transform spectrom- 
eter (EFTS) backends jRisacher et al.ll2006t iKleine t al. 2006j). 
The signal and image sidebands are separated by 12 GHz. We 
used two setups in our observations, which are described in 
Tabled Each EFTS has a bandwidth of 1 GHz, and 8192 
channels, which corresponds to a velocity resolution of ~0.12- 
0.15km s-i for lines between 285-350GHz. For the H2CO(4- 
3) rotational transition, we observed only the 4o4 - 3o3 K level, 
since the remaining K levels of this transition were outside the 
bandpass due to a mistake in the center frequency. 

Each observation was performed with on-off iterations us- 
ing three subscans, integrating in total between 90-230 seconds 
on source. From the IRAM 30m observations, which were car- 
ried out before the start of our APEX runs, it was clear that the 
off positions were often contaminated by extended CO emis- 
sion. Therefore, we used for the APEX observations off posi- 
tions much further out, at 1800", which in most cases were free 
of emission. 

The calibration of the APEX data was, just as with the IRAM 
30m telescope, carried out using the chopper-wheel technique. 
The telescope efficiencies, necessary to calculate the main beam 
brightness temperature, are listed in Tablefefta:freq. Before each 
run, the focus was adjusted on a planet, usually Jupiter Pointings 
were made once every three hours on nearby Hii regions with 
strong dust continuum emission, such as SgrB2(N). The precip- 
itable water vapor was between 1 .44 and 2.80 mm, while the sys- 
tem temperatures ranged from 190-300 K at 280 GHz and 410- 
575 K for 356 GHz. 

2.4. Data reduction 

The processing of the data, such as smoothing (spectral averag- 
ing) to increase the signa l-to-noise r atio and baseline subtrac- 
tion, was done in CLASS (lPetvll2005h . For most lines, Gaussian 
fitting was performed to retrieve the basic line parameters, such 
as line widths, vlsr, and line intensities. For the lines with hy- 
perfine structure (hfs), we used the hfs method, which allows 
additionally the derivation of the optical depth of the main hfs 
component (when applicable). The spectral plots were also pre- 
pared in CLASS. 

3. Results 

An overview of all the detected lines per clump and the clump 
J2000 positions are given in Table IB. fl For sources with no de- 
tection, we put a 3cr upper limit (also in Table IB. it . 

3.1. Infali 

The infall signature of a source can be recognized by a line pro- 
file with a double-peaked structure, where the intensi t y of the 
blue peak exceeds the intensity of the red peak dLeung & BrowrJ 
1977tlZhou et al.ll993tlMvers et alJl996HMardones et al.ll997l: 
Evan sl ll999l) . Infall models assume a source where the infall ve- 
locity, density, and excitation temperature increase toward the 
center. In this scenario, the blue peak of an optically thick self- 
absorbed line originates in the rear part, and the red peak in the 
front part of the infalling shell. An increase in the infall velocity. 
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Table 1. Properties of clumps in high extinction clouds based on Paper I 



diffuse clouds peaked clouds multiply peaked clouds all clouds 







mean 


range 


mean 


range 


mean 


range 


mean 


cloud mass 


Mo 


495 


17-3039 


1420 


70 - 5500 


1900 


500 - 6500 


910 


cloud diameter 


pc 


0.8 


0.2-2.3 


1.7 


0.7 - 3.7 


1.9 


1.2-3.8 


1.2 


clump mass 


Mo 


105 


12 - 283 


130 


28-738 


185 


21-431 


150 


clump A'ht 


10^2 cm-2 


5.4 


2.9-9.5 


5.5 


3.3-7.7 


8.5 


2.5-26.3 


6.8 


clump diameter 


pc 


0.3 


0.11-0.56 


0.33 


0.14-0.76 


0.36 


0.15-0.47 


0.3 


clump Av 


km s"' 


1.2 


0.7-2.3 


1.4 


0.9-2.5 


1.6 


0.8-2.8 


1.5 


clump d 


kpc 


2.8 


1.1-5.7 


3.4 


1.9-7.2 


3.6 


2.1-4.7 


3 


clump temperature 


K 


13.5 


9.3-16.7 


15.7 


11.9-18.6 


17.5 


12.4 - 24.7 


16 


clump E 


gcm-- 


0.3 


0.2-0.5 


0.3 


0.2 - 0.5 


0.5 


0.2-1.5 


0.4 



Notes. Rows are (from top to bottom) cloud mass, cloud diameter, clump mass, clump column density, clump diameter, clump line width, clump 
kinematic distance, clump rotational NH3 temperature, and clump surface density. 



Table 2. Molecular lines and frequencies 



Molecule 


J+\^ J 


EJks 


Frequency 


ncri(° Bandwidth 


Resolution 


Beam 




-fcff 


Icr r.m.s 






(K) 


(GHz) 


(cm-3) 


(MHz) 


(km s-') 


(") 






(K) 








IRAM 30m telescope 














H"CO+ 


1-0 


4.16 


86.754 


1.7 X 10' 


160 


0.16 


29 


0.78 


0.98 


0.06 


SiO 


2-1 


6.25 


86.847 


7.3 X 10' 


160 


0.16 


29 


0.78 


0.98 


0.07 


HCO+ 


1-0 


4.28 


89.189 


1.8 X 10' 


40 


0.08 


28 


0.78 


0.98 


0.10 


CH3CN 


5^, K=Q, 1, 2, 3 


13.24* 


91.987* 


4.7 X 10' 


120 


0.16 


27 


0.78 


0.98 


0.09 


N2H+ 


1-0 


4.47 


93.174 


1.6 X 10' 


40 


0.08 


27 


0.78 


0.98 


0.11 


C'^'O 


2-1 


15.81 


219.560 


9.2 X 10^ 


512 


1.5 


12 


0.62 


0.94 


0.12 


CO 


2-1 


16.60 


230.538 


1.1 X lO'* 


512 


1.5 


11 


0.58 


0.92 


1.3 








APEX telescope 














N2H+ 


3-2 


26.83 


279.512 


3.0 X 10* 


1000 


0.15 


22 


0.73 


0.97 


0.23 


H2C0 


404-303^^ 


34.90 


290.623 


9.2 X 10* 


1000 


0.15 


22 


0.73 


0.97 


0.20 


HCO+ 


4-3 


42.80 


356.734 


9.1 X 10* 


1000 


0.12 


18 


0.73 


0.97 


0.47 


CO 


3-2 


33.19 


345.796 


3.5 X 10* 


1000 


0.12 


18 


0.73 


0.97 


1.1 



Notes. Columns are (from left to right) the molecule, its transition, its upper energy level, its frequency, its critical column density, the bandwidth 
used in the observation, the velocity resolution of the observation, the telescope beam, the main beam and forward beam efficiency of the obser- 
vation , and the mean Icr noise value. Calculated from the collision rates at T = 20K from the LAMBDA molecular database dSchoier et al.l 
120051) . Frequency and E^/kg for the K=Q level. Observed were ^"=0, 1, 2, 3, 4 The observed K level is given in subscript. 



Vin, can cause the red peak to diminish; at very large values of Vin, 
the red peak can even disappear and become a red "shoulder". A 
red excess, in contrast to this blue excess, may be caused by ex- 
pansion or outflow. Alternatively, the red excess can be caused 
by an outwards moving blob of matter, instead of large-scale 
outward motion (lEvanslll999h . To interpret the infall profile of 
an optically thick line, it is also necessary to compare it with 
an optically thin line measured toward the same position. The 
latter, often a line from a rare isotopologue, shows a maximum, 
defining the systemic velocity at the self-absorption minimum 
of the optically thick line. The observed line profiles can, in ad- 
dition to large scale motions such as infall, also be influenced 
by abundance changes though the clouds. Hence, it is desired 
to investigate infall using several molecules. To probe different 
depths in the clouds and thereby possible different infall veloci- 
ties, we used lines with different critical densities: the CO(2-l), 
CO(3-2), HCO+(1-0), and HCO+(4-3) transitions. 



3.1.1. CO line profiles 

We started the study with the CO(2-l) and CO(3-2) lines, whose 
critical densities are around IxlO"* and 4x10"* cm""* (Table|2]l, re- 
spectively. Since the CO molecule, as a common component of 
the interstellar medium (ISM); has a very extended distribution, 
the off position is not necessarily free of its emission. This can 
result in artificial absorption lines in the spectrum. To avoid this, 
one can observe with two different off positions and/or use dif- 
ferent CO transitions. Since higher / transitions have a higher 
critical density, their emission is confined to a more compact 
and higher density region than emission from lower J lines. We 
employed this strategy and observed, in addition to the CO(2-l) 
line with the IRAM 30m, the CO(3-2) line with APEX using dif- 
ferent off positions. Also, because of the molecules' ubiquity and 
relatively low critical density, CO observations pick up emission 
from various Galactic arms. This manifests itself in the spectrum 
by more than one peak at various local standard of rest (LSR) ve- 
locities. 

We searched for infall signatures by comparing the veloci- 
ties of the CO(2-l) and CO(3-2) emission peaks with those of 
the optically thin J8o(2-l) line. When the peak of the optically 
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; H CO (1-0) 

HCO* (1-0) 



50 

Velocity (km/s) 

Fig. 1. Two examples of infall/outflow signatures of the '^CO 
line profiles. The optically thin J8o(2-l) emission (scaled for 
better visibility, shown in green) indicates the systemic velocity. 
The CO(2-l) and CO(3-2) emission is shown in black and red, 
respectively. The panels show infall (top ) and outflow (bottom). 



thick line of the main isotopologue appeared shifted blueward of 
the J8o peak, the source was marked a blue excess source or a 
blue source. A source with a red shifted peak was marked a red 
source. 

In general, both the CO transitions showed similar behavior 
(see Table IBT2I 1. Two examples of infall and outflow signatures 
are presented in Fig.lT] When the peaks of the optically thick and 
thin lines were at the same velocity, nothing could be inferred. 
Mostly, the CO lines were broader than the J 80 lines and had 
line wings extending from 10 to almost 50km s"' relative to the 
systemic LSR velocity. Line width broadening is partially due to 
a high optical depth, which is higher for the CO lines than for 
J80 lines. However, the wings are dominated by gas undergoing 
strong and dominant motions, such as infall or outflow, and in 
cases such as G022.06+00.21 MMl by high-velocity emission 
features ('bullets'). 

The results for CO are given in Table lB.2l For the lower tran- 
sition, the presence of wings is indicated, and, if present, their 
velocity range is given. The CO line profiles delivered the first 
indications of infall or outflow. A deeper study was done with the 
HCO+(1-0) and H'''CO+(1-0) fines, reported in the next section. 

3.1.2. The skewness parameter 5v 

The optically thick HCO^(l-O) line profiles were used for an 
in-depth study of infall and outflow motions in the clumps. The 
systemic velocity of a clump was determined from the optically 
thin Hi3CO+(l-0) fine. Examples of HCO+(1-0) fine profiles 
are shown in Fig. |2] For half of the clumps, the HCO^(l-O) 
line showed a double-peaked profile with a self-absorption dip. 
Here, the ratio of the two intensity peaks, ^p^, was used to check 

■' red 

whether the source showed blue or red excess in the emission 
(ratio listed in Table lB.2l . 

When the HCO^(l-O) line profile showed only one peak, 
we compared its velocity with the velocity of the optically thin 
Hi^CO+(l-0) peak. The measure of this shift in velocity is 
ca lled the skewne s s par ameter Following the method outlined 
in lMardones et al.l (119971) . the skewness parameter 6v, is defined 



GOI 4.63-00.57 MMl 




50 55 
Velocity (km/s) 

Fig. 2. Three examples of infall/outflow signatures of the 
HCO+(1-0) line profile. Red marks both Hi3CO+(l-0), which 
is optically thin, and the systemic velocity of the dense gas. 
Compared to this the shift of the peak of the HCO^(l-O) emis- 
sion, in black, becomes clear The panels show infall (middle 
and bottom), outflow (bottom right), and a case of central self- 
absorption, where both peaks are equal (top). The profiles show 
various degrees of self-absorption. 



as 



Vthick - Vthh 



Av 



(2) 



thin 



where Vthick and Vthin are the LSR velocities of the peaks of the 
optically thick and thin lines, respectively; Avthin is the line width 
of the optically thin line. The line width and LSR velocity of 
the H'^CO^(I-O) line were retrieved by Gaussian fits (Table 
IB3T i. For the optically thick HCO+(1-0) line, the profiles of- 
ten showed non-Gaussian shapes; therefore, the position of the 
peak was determined by eye (also given in Table IB.3l l. After 
Ma rdones et al.l (Il997h . we adopt the definition of significant 
blue and red excess, namely for the former if 6v < -0.25 and 
the latter when 6v > 0.25. 

The skewness parameter was determined for the HCO^(l- 
0), HCO+(4-3), and CO(3-2) fines (afi fisted in Table iR^l i. 
For all molecules, the values ranged between -1.5 and 1.5, ex- 
cept for G14.39-00.75B where the Sv of the HCO+(1-0) fine 
was -5.79. The latter is due to a very blue shifted HCO^(l-O) 
peak, most likely caused by a high self-absorption in this source. 
The clumps in high extinction clouds have skewness parame- 
ters similar to values i n the literature: -0.5 for U CHii regions 
( Wv rowski etaP 120061 -1.5 to 1 for HMPOs (Fuller etaL 
20051), -1 to 2 for Extended Green Objects (EGOs, Chen et^L 
2010), which are sho cked regions, a nd -0.5 for the cluster- 
forming clump G24.4 (IWu et al.ll20"05l) . 
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Table 3. Distribution of the skewness parameter per molecule 

3 - 



Transition 




A'.-ed 


-^skew 


A'.o, 


E 


P 


< <5v > 


HCO+(1-0) 


24 


14 


38 


47 


0.22 


0.07 


-0.26 


HCO+(4-3) 


5 


4 


9 


18 


0.06 


0.50 


-0.11 


00(2-1)" 


24 


11 


35 


35 


0.37 


0.02 




CO(3-2) 


8 


7 


15 


16 


0.06 


0.50 


-0.31 



Notes. Columns are (from left to right) number of blue excess sources, 
number of red excess sources, total number of sources with a significant 
skewness, total number of sources, the excess parameter, the probability 
of the distribution to arise by chance, and the average skewness. 

The numbers of the C0(2-l) transition serve only as an indication. 

Cross correlations of the skewness between the three differ- 
ent molecules (Fig.|3] clumps from peaked and multiply peaked 
clouds are color coded in light of the evolutionary sequence dis- 
cussion in Sect. 14.2b show that the HCO^(4-3) line is the least 
sensitive to infall and/or outflow, while the HCO^(l-O) seems 
the most sensitive. Most of the HCO^(4-3) lines' 6v values are 
< 10.251, while those of the HCC >+(l-0) an d CO(3-2) lines ex- 
hibit a range of skewness. Also. iFuller et a l. (2005) found that 
the transitions with higher critical densities show less infall. 
These observations suggest that the infall occurs predominantly 
in the low-density environment, which is in contrast with the 
increasing infall velocity toward the center that is expected in 
the collapsing clump model. The observations can possibly be 
explained by the different optical depths of the HCO^ transi- 
tions, because an infall profile requires the transition to be suf- 
ficiently opaque. Depending on the density of the infall envi- 
ronment, the transitions with a high critical density might not be 
opaque enough to observe the skewness profile, while transitions 
of lower critical density will be already sufficiently optically 
thick to observe the infaU profile. For the HCO+(1-0) and CO(3- 
2) lines, the skewness seems to be correlated, which makes the 
infall determination based on these two molecules more reliable. 

The excess parameter shows the average behavior of the blue 
excess over the red excess sources. It is defined as £ = (A^biue - 
A^red)/Mot. where A^i-ed is the number of red excess sources, A^biue 
the total numb er of blue excess sourc es, and A^tot the total number 
of all sources (iMardones et al.|[T997h . 

Table |3]lists the excess parameter and average skewness pa- 
rameter by transition. For each transition, we calculated, through 
the binomial test, the probability P that the distribution between 
blue and red excess sources occurred by chance. The binomial 
distribution is defined as 

p = i^^Ai-py"-''\ (3) 

where n is the total number of trials, k the number of successes, 
and p the success probability. In our case, n is the total number of 
blue and red excess sources, k the number of blue excess sources, 
and the success probability p equals 0.5 if the blue and red ex- 
cess sources are randomly distributed. We can then calculate the 
possibility that the distribution of the number of blue sources 
equal or higher than the number observed arises by chance by 
adding all possibilities P(n, k, p) + P{n, k+1, p) + .. until k - n. 

The small probability P of 7% for HCO+(1-0) fine in Table 
[3] indicates that there are significantly more blue excess sources 
than expected for a random distribution. On the other hand, the 
higher CO(3-2) and HCO^(4-3) transitions have a probability 
of 50%. Hence, the number of blue excess sources is likely by 
chance and not significant. The results for the CO(2-l) line. 
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Fig. 4. N2H+(l-0) {top) and N2H+(3-2) (bottom) line profiles. 
The green lines show the locations of the hyperfine components. 



P - 0.02, should be treated with care, because the same defi- 
nitions of blue and red sources do not apply. Due to the possible 
confusion in the line profiles, the CO(2-l) was just classified 
based on the peak emission with respect to the J8o(2-l) line 
without the calculation of Sv. Therefore, the numbers in Tables 
IB.2l andl3]have to be taken just as indications. In conclusion, the 
HCO^(l-O) line shows a significant blue excess and is the best 
indicator of infall for our sample of clumps. The excess param- 
eter of the HCO^(l-O ) line, 0.22, is s imilar to values found in 
a previous survey by F uller et aP (12005) of high-mass protostel- 
lar objects reaching excesses of 0.29 and 0.31. Studies of low- 
mass star formation also report similar numbers (IMardones et al.l 
1997: Evans..2003.) . 

3.2. Temperature and density estimations from N2H^(1-0) 
and (3-2) transitions 

The N2H+(l-0) and N2H+(3-2) rotational transitions have hfs 
arising from the interaction between the molecular electric field 
gradient a nd the ele ctric quadrupole moments of the two nitro- 
gen nuclei dCaselU et al. 1995). The N2H+(l-0) rotational transi- 
tion is split into seven hyperfine components. The three main hy- 
perfine groups were well resolved in our observations (see Fig. 

and for a few sources we could fit all seven components. In 
the higher transition, we were barely able to even resolve the 
three main groups of hfs components (Fig.|4|. 

For fitting the N2H^(l-0) transition we took the hfs into ac- 
count by using the hfs method for the line fitting in CLASS B 
This method assumes one excitation temperature and one line 
width for all seven hyperfine components, as well as the fact that 
the opacity as a function of frequency for each hfs has a gaus- 
sian shape. Then, besides the intrinsic line width and integrated 
intensities, the hfs method also determines the optical depth of 
the main hfs group. The total optical depth Ttot can be calculated 
from the obtained main hfs group optical depth Tmain by mul- 
tiplying it by the sum of relative intensities of the satellites S 
as Ttot = '^Tmain- In casc the relative intensities are scaled so that 
their sum equals unity, the hfs fit directly delivers the total optical 
depth (the sum of all the hfs optical depths). We calculated the 



" http://www.iram.es/IRAMES/otherDocuments/postscripts/classHFS.ps 
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(5v HCO^(l-O) 6v HC0^(4-3) 

Fig. 3. Skewness parameter of the HCO^(l-O) line versus that of the HCO^(4-3) and '^CO(3-2) lines. The dashed lines mark the 
boundary of significant excess at 6v > |0.25|. The blue triangles represent clumps in peaked clouds, while the red squares represent 
sources in multiply peaked clouds. 



colum n density of N2H^, A^N2H+ ( IA.3I ) following 'Bens on et alJ 
(119981) . The results of the N2H+(l-0) Hne profile fitting; the inte- 
grated intensity (including the hyperfine components), total op- 
tical depth, and line width are placed together with the excitation 
temperature and the column density in Table IB. 4l For N2H^(3- 
2), Table IB. 5] lists the parameters obtained by Gaussian fits. We 
did not observe the hfs as clearly as in the lower transition and 
could therefore not determine the optical depth for the higher 
transition. 

For sources with r(i_o),tot ^ 0.4, which usually also had high 
relative errors, we found unrealistically high Tex values that were 
much higher than the ammonia rotational temperatures, T^ot 
(Paper I). The CLASS hfs method does not provide good esti- 
mates for excitation temperatures (hence also column densities) 
when the relative errors in optical depth are very large (> 40%) 
and when the optical depth is very small (t(i_o) tot ^ 0.4) (see, 
e.g., iFontani et alJl2006l l201 Ih . After discarding these sources, 
the excitation temperatures generally ranged from 3 to 29 K, av- 
eraging 8.1 K . In fact, only two sources had a Tex > 20 K; 
G017. 19+00.81 MM2 (28 K) and MM3 (29 K). We assumed 
that the N2H^ gas has a similar kinetic temperature as the 
NH3 gas and that they are both in local thermal equilibrium 
(LTE), in which case the Tgx should equal Trot. However, the 
determined excitation temperatures were, with the exception of 
G017.19-h00.81 MM2 and MM3, all lower than the rotational 
temperature from ammonia. The diff'erence between the Tex de- 
termined here and the Tiot determined by ammonia is the N2H^ 
filling factor While Tgx, is inversely proportional to the filling 
factor (see lA.ll ). the rotational temperature is an intrinsic tem- 
perature (independent of / since it is solely determined by line 
ratios). We estimated the filling factor from the ratio of the ex- 



citation temperature of N2H^(l-0) and the NH3 rotational tem- 
perature: 

7'ex,N2H+ 

/N,H+(l-()) ~ — ■ (4) 

J rot, NH, 

The individual filling factors range from 0.2 to 0.8. For clumps 
G17. 19+00.81 MM2, and MM3, we found meaningless filling 
factors larger than unity, which is possibly due to an overes- 
timation of the excitation temperature. A filling factor smaller 
than unity indicates either that the source size is smaller than the 
beam or that the source is clumpy instead of centrally condensed. 
Alternatively, the N2H^(l-0) line could be sub-thermally ex- 
cited, a consequence of densities lower than the critical density. 
At such densities, a transition's upper energy level becomes un- 
derpopulated: collisions do not manage to populate it, and the 
system is in non-LTE. In this case, the excitation temperature 
will be lower than the kinetic temperature, thus mimicking the 
behavior of a small filling factor. 

To have another estimate of the filling factor, we performed 
radiative transfer modeling to obtain the intrinsic N2H^ column 
density. In a non-LTE system, the ratios of line intensities devi- 
ate from the LTE predictions and become a function of density, 
column density, temperature, and optical depth. In particular, 
we consider the N2H^(l-0)/N2H^(3-2) line ratio and combine 
t his with non-LTE rad iative transfer modeling using RADEX 
(van der Tak et al. 2007) and the molecular database LAMBDA 
(Schoier et al..2005) . Figure|5]shows the integrated intensities of 
the two N2H"^ transitions plotted against each other We use dif- 
ferently colored symbols for each cloud class (diffuse, peaked, 
and multiply peaked) in the graphs shown throughout this paper 
to visualize the (possible) behavior for different phases of cloud 
evolution. The integrated intensity, j Tmb'^v, of the lower tran- 



7 



K. L. J. Rygl et al.: Infall and onset of star formation 



sition, N2H^(l-0), exceeded that of N2H^(3-2) by an average 

JWv(N2H^(l-0)) ^ 

/ rMBdv(N2H+(3 - 2)) 

In this ratio, we assumed that the filHng factor ratio is 1, since it 
is hkely that the filHng factors of the two transitions cancel out 
as the APEX beam is only slight smaller than the IRAM 30m 
beam (Table |2] lists all beamsizes). There are three sources that 
lie far from the average: G014.63-00.57 MMl, G017.19+00.81 
MM2, and G022.06+00.21 MMl, which are all bright at 24;um 
and contain water masers (see Table IB. lb . 



The RADEX results for the volume density were almost all 
within a factor two or less of the measured volume density from 
the 1.2 mm continuum (Table |4]i. Given the uncertainties in the 
derivations of both volume densities, we can say that they are 
in agreement for most of the sources. The N2H^ ratio and the 
N2H^(l-0) optical depth can therefore be used to estimate the 
hydrogen volume density and intrinsic N2H^ column density of 
a clump for a given kinetic temperature. 

We estimated the abundance, of N2H^ by comparing the 
NnjU* corrected by the filling factor from the temperature ratio 
(for each source individually) to the hydrogen column density, 
A^H2 > measured from the 1 .2 mm continuum (Paper I): 



50 



-10 



30 



20 



10 - 




(6) 



1 20 
/TmbAv N2H^(3-2) [K km 



30 



Fig. 5. N2H+(l-0) integrated intensity versus N2H+(3-2) inte- 
grated intensity for all detected sources. The blue triangles rep- 
resent clumps in peaked clouds, while the red squares represent 
sources in multiply peaked clouds. The green solid line marks 
the trend of the observed ratio of 2.9. 

The observed N2H^ column density is roughly proportional 
(depending on the value of Tex, see Eq. IA.3l l with excitation 
temperature, A^n,h+ for 5K< Tgx <10K, which in turn 

is inversely proportional with filling factor Hence, the ratio of 
observed N2H^ column density (from the N2H^(l-0) line) and 
the intrinsic column density calculated by RADEX (based on 
the N2H+(l-0)/N2H+(3-2) integrated intensity ratio) can be in- 
terpreted as a result of the filling factor of the N2H^(l-0) line. 
Using RADEX, we estimated the behavior of A^n,h+/Av and «h, 
for different values of the N2H+(l-0)/N2H+(3-2) Hne ratio and 
the opacity of N2H^(l-0) for kinetic temperatures of 10, 15, 20, 
25 K. We chose this temperature range based on rotational tem- 
peratures derived from the ammonia observations (see also Table 
[T]). The results are shown in Fig.|6]and TableH) 

In Fig.|6]we marked the positions of the sources from Table|4] 
at the intersection of the ratio of the (1-0) and (3-2) N2H^ inten- 
sities and the N2H^(l-0) optical depth in the plot calculated with 
a kinetic temperature closest to the NH3 rotational temperature 
of the clump. To obtain the intrinsic N2H^ column density, we 
multiplied jVnjH+ZAv by the observed N2H+(l-0) line width. In 
the plots, we see that a rising Tidn increases the column density 
for a given optical depth and a line ratio, while the volume den- 
sity decreases. For the handful of sources that we have measured 
both N2H^ transitions, we compared the filling factor based on 
the column density with the one determined from the temper- 
ature. On average, the filling factors agreed within 30%. Here 
we excluded clumps G017. 19-1-00.81 MM2 and MM3, since col- 
umn density-derived filling factor came out larger than unity, and 
which confirms that for these sources the derivation of the N2H^ 
parameters, from which the excitation temperature and column 
density are calculated, gave unrealistic values. 



The individual clump abundances are listed in Table IB. 41 and 
range from 0.26 - 4.0 x lO"'', with a mean of 1.3 x 10 . 
The true abundances are likely to be slightly lower, since the 
A^H2 is a beam-averaged column density and contains an un- 
known filling factor. Taking this into account, t he ab undance is 
in agreement with the value of iPirogov et al.l (|2007|) for high - 
mass star-forming cores and the results of Ragan et aP (l2006l) : 
IVasyunina et al.. (2011.) for IRDCs. 



3.3. Presence of young stellar objects 

3.3.1 . Hot molecular cores: CH3CN and H2CO 

CH3CN and H2 CO are usually found to\yard YSOs or hot 
molecular cores jMangum& Woottenlll993t lObni et al.lll996t 
Ivan der Tak et al. 2000*) , which are wa rm (> 100 K) and dense 
(> 10^ - 10''cm-3) (K urtz et al.ll2000h . Both molecules can be 
used as a tracer of ongoing star formation. 

CH3CN is a symmetric-top molecule where each rotational 
level is split for different projections of angular momentum 
along the symmetry axis of the molecule labeled by the quantum 
number K. Within one rotational level, 7-1-1 — > 7, the transitions 
of the K components are determined solely by collisions, and 
their relative intensities are therefore relate d to the kinetic tem- 
perature of the region dSolomon et aLlfl^TTb . The K components 
with one rotational level are very closely spaced in frequency 
and can be observed simultaneously with one backend with a 
bandwidth >30MHz. Under the assumption that the CH3CN(5- 
4) A'=0, 1, 2, 3, and 4 emission is optically thin and uniformly 
fills the beam and that all K transitions within this J level are 
characterized by one temperature, we derived the rotational tem- 
peratures, T mt, and the total CH3CN column density following 
lArava etak l (12005). These authors studied the CH3CN emission 
toward objects similar to ours (massive star-forming clumps) and 
showed that it is reasonable to assume the optically thin limit for 
this line. A beam-filling factor smaller than unity would affect 
the K levels similarly and hence would have no large influence 
on the temperature derivation if the filling factor does not vary 
from one K component to another. The CH3CN column den- 
sity, however, is affected by the beam-filling factor Given the 
beam size of 27" against the average angular clump size found at 
1.3 mm, ~20", the beam-filling factor could be slightly smaller 
than unity. 

We detected the K=Q to 3 levels, or a subset of these, for 
several sources; the data were too noisy (1 cr rms of 0.9 K) for 
a K-4 level detection. All detections, derived r,ot, and column 
densities are given in Table IB. 61 The CH3CN rotation temper- 
atures were much higher than the NH3 temperatures (see Table 
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Table 4. Results of RADEX calculations based on the (l-0)/(3-2) N2H+ ratio 



N2H+(l-0) N2H+(3-2) 



Source name 




1 T^AJudv 

J ^ MB"^ 
t^JV Kill S ) 


Av 

^^KIIl S ) 


7" 


1 T\/id(\v 
J ^ MB*-*'^ 

(^JV Kill S ) 


^ rot 


cm ) 


N 

cm ) 


''Nht .obs 

cm ) 


'•Nhj ,r 

\l\J CIIl ) 


G013.91-00.51 


MMl 


12.9 


1.19 


5.34 


2.4 


14.1 


2.2 


5.8 


2.2 


0.3 


G014.63-00.57 


MMl 


45.8 


2.46 


2.35 


25.6 


18.1 


8.9 


18.7 


7.7 


3.4 


GO 14.63-00.57 


MM2 


17.9 


1.27 


3.90 


8.2 


15.7 


2.8 


9.4 


7.3 


3.4 


G017. 19+00.81 


MM2 


32.9 


2.21 


1.0 


18.7 


18.7 


13.5 


7.7 


6.3 


4.0 


G017. 19+00.81 


MM3 


25.1 


1.74 


0.91 


1.7 


20.1 


10.9 


2.7 


0.8 


1.8 


G018. 26-00.24 


MMl 


23.5 


2.50 


2.10 


8.6 


18.2 


3.7 


10.5 


1.8 


1.4 


G018.26-00.24 


MM2 


26.0 


2.16 


4.00 


7.6 


17.4 


4.3 


10.1 


1.3 


1.0 


G022.06+00.21 


MMl 


28.7 


2.50 


1.17 


16.3 


24.7 


7.7 


13.0 


6.8 


5.2 


G024. 37-00. 15 


MM2 


12.7 


2.28 


1.55 


5.4 


15.7 


2.1 


6.2 


1.2 


2.9 


G024.94-00.15 


MMl 


20.0 


2.48 


1.30 


7.4 


15.2 


4.2 


5.0 


2.0 


2.9 


G030.90+00.00A 


MMl 


22.8 


2.38 


0.26 


1.1 


18.6 




3.1 


0.7 


1.8 


G034.7 1-00.63 


MMl 


11.7 


2.31 


2.61 


6.1 


17.8 


1.9 


17.6 


1.0 


3.4 


G034.7 1-00.63 


MM2 


11.5 


1.58 


4.90 


2.7 


12.4 


2.1 


6.3 


1.1 


1.4 


G053. 81-00.00 


MMl 


9.0 


1.64 


0.56 


3.7 


12.4 


1.0 


12.1 


3.6 


8.6 



Notes. Column are (from left to right) source name, integrated N2H^(l-0) line intensity, FWHP N2H+(l-0) line width, N2H+(l-0) optical depth, 
integrated N2H^(3-2) line intensity, NH3 rotational temperature (Paper I), observed N2H^ column density, calculated N2H^ column density from 
RADEX, observed hydrogen density from the 1.2 mm continuum (Paper I), calculated hydrogen density from RADEX. 



umn density of 5 x 10^^ cm ^, almost all clumps show infrared 
emission and signs of outflows. 

For both the CH3CN and H2CO molecules, there is a gen- 
eral cut-off in the column density at 4 x 10^^ cm"^, below 
which no YSOs are found. This cut-off is slightly lower than 
the column density of ~1 x lO^^ cm'^ theoretically require d 
to form massive stars (0.7 gem -2, iKrumholz & McKeell2008l) . 
iLopez-Sepulcre et al] (|2010) found a colunm density threshold 
of 7.7 X 10^2 cm"2, below which the outflow rate decreases 
rapidly and the outflows are less massive. Hence, while obser- 
vational data seem to strongly support a column density cut- 
off for (massive) star formation, they also seem to indicate a 
slightly lower threshold than theoretically advocated. The dif- 
ference, however, is not too large, and the observed column den- 
sities are beamaverages, which means that if the source size is 
smaller, we are underestimating the true column density. 

3.3.2. Outflow: SiO(2-1) emission 

Apart from evidently heating their envi ronments, (massive ) 
YSOs drive strong molecular outflows (iBeuther et al.l |2002|) . 
which may remove, in the case of disk accretion, the angular 
momentum away from the forming star. While mechanisms to 
accelerate the outflow have been proposed for low-mass star for- 
mation (Shu et al. 2000), in the high-mass scenario the accelera- 
tion mechanism is not well understood yet. 

Interstellar dust grains contain silicates (lDrainel2003h . Shock 
waves driven by outflows can sputter Si-bearing compounds in 
dust grain s leading to increas ed SiO abundances in the post- 
shock gas (ISchilke et al.lll997l) . Thus, SiO is expected to be as- 
sociated with ongoing star formation via its formation history. 
Observations support this, since SiO is more common in more 
active (or evolv ed) sources, while very wea k in quiescent (hence 
early) sources dBeuther & Sridharanll2007h . Naturally, one also 
expects to observe wide line widths. We found line widths from 
~ 1.5-7 km s"' for very weak SiO detections, while for stronger 
SiO emitting sources the line widths increase up to 26 km s"'. 
The individual values for all detections are given in Table IB. 71 
along with the main beam brightness temperatures and LSR ve- 
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Fig. 7. CH3CN rotational temperatures versus NH3 rotational 
temperatures. Clumps for which the CH3CN temperature had er- 
rors smaller than 12 K are red. 



[]and Fig.|7]), since the CH3CN traces warmer gas than ammo- 
nia does. For the H2CO(4o3-3o4) line, the observed main beam 
brightness temperatures, LSR velocities, and line widths were 
determined from Gaussian fits and are listed in Table IB. 7l 

To understand if YSOs are preferably present in denser or 
more compact clumps, we plotted all the clumps with evidence 
for an embedded YSO in a size versus the column density graph 
(Fig. [8]). In addition to plotting clumps with evidence of a hot 
molecular core through the detection of H2CO and CH3CN, we 
plot clumps with a 24 /im detection (direct evidence of dust 
heated by an YSO) and SiO emission (evidence of outflows, 
see Sect. I3.3.2l i. Apparently, YSOs are found in clumps of all 
sizes, but, more interestingly, they are mostly found at high 
column densities. The minimal A^Hj at which the H2CO(403- 
3o4) fine was detected is 4 x lO^^cm"^, while the CH3CN(5-4) 
K - Oemission was still observed toward one less dense clump 
with A^H, = 3 X 10^2 cm"2. The column density cut-off is most 
clearly seen in the 24/j.m and SiO line detections; above a col- 
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Fig. 8. High extinction clumps placed in a diagram of column 
density versus the FWHM diameter (in pc) derived from the 
1.2 mm continuum (Paper 1). The detections of MIPS 24fim 
emission, H2CO(4o3-3o4), CH3CN(5-4) K^Q, and SiO(2-l) are 
marked in color, black symbols are non-detections. The crosses 
mark the clumps in diffuse clouds, triangles clumps in peaked 
clouds, and squares clumps in multiply peaked clouds. 
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Fig. 9. SiO detections against the HCO^(l-O) skewness param- 
eter 6v. The green crosses represent sources in diffuse clouds, 
the blue triangles clumps in peaked clouds, and the red squares 
sources in multiply peaked clouds. 



locities. The mean noise level of the SiO observations is 0.07 K, 
which is comparable to previous studies in massive star-formin g 
regions by Beuther & Sridharan [200% and lMotte et all (l2007h . 
which reached noise levels of 0.02 K and 0. 10 K, respectively. 

We detected the SiO(2-l) transition toward 50% of the 
whole sample. Of all the SiO-emitting sources, a remarkable ma- 
jority, 70%, coincided with infall sources (having a HCO^(l-O) 
skewness parameter, 6v < -0.25), as illustrated in Fig. |9] This 
finding indicates that accretion is accompanied by collimated 
outflows, which cause the enhancement of SiO. A further 15% of 
the SiO detections are sources exhibiting expansion (6v > 0.25), 
and the remaining 15% are sources showing neither expansion 
nor infall (-0.25 < Sv < 0.25). Clumps without SiO detections 
have HCO^(l-O) skewness parameters divided almost equally 
between infall (ten clumps) and expansion (eight clumps). 




Fig. 10. SiO line width against the hydrogen column density 
A^H2- The green crosses represent sources in diffiise clouds, the 
blue triangles clumps in peaked clouds, and the red squares 
sources in multiply peaked clouds 



To check if the SiO line widths, and hence the outflows, 
are stronger for the more evolved sources, we compared them 
with the hydrogen column density determined from the 1 .2 mm 
continuum. Figure [TO] shows that up to column densities of 
10^^ cm"^ there is a tendency of increasing line widths with 
column density. The few clumps that are above this column 
density limit, G014.63-00.57 MMl, G017.19-h00.81 MM2, 
G022.06+00.21 MMl, and G014.63-00.57MM2, show no clear 
trend. The SiO emission was detected toward all multiply peaked 
clouds, but just toward three peaked clouds (25%) and only two 
diffuse clouds (18%), singling out the multiply peaked clouds as 
active regions of star formation. 



3.3.3. Prestellar and starless sources 

Our sample contained 19 sources with no evidence of a YSO. 
Of these sources, one-third shows infall in the HCO^(l-O) line 
profile indicating that these are in a prestellar phase, one-third 
is expanding according to the HCO^(l-O) skewness parameter, 
and the last third shows neither infall nor outflow. The sources 
wit hout any sig n s of in fall or YSOs were defined as starless. 

iMotte et alJ (l2007l) showed how difficult it is to find high- 
mass analogies of lo w-mass starless cores, such as L1521F 
(IShinnaga et al.ll2004l) . in their study of the Cygnus X complex. 
Remarkably, they found SiO emission in all of the observed 
infrared-quiet massive dense cores, indicating that star forma- 
tion is already ongoing. 

In our sample, we found six prestellar and six starless 
sources. Only half of them contain a clump (an overdensity of 
typical hydrogen column density of 10^^ cm"^, which could be 
identified by the Miriad routine "sfind", see Paper 1); most of 
them are too diffuse in the mm continuum. The masses of the 
prestellar clumps, derived from the 1.2 mm continuum, were 
greater than their virial masses, indicating that gravity is the 
dominant force. However, also for the starless clumps, their 
virial parameters indicate that gravity dominates their internal 
dynamics. Apparently, it is not merely the lack of gravitational 
force that separates the starless from the prestellar clumps. The 
prestellar and starless sources were found in the diffuse and 
peaked clouds, but not in the multiply peaked clouds, which 
agrees with the previous indications of the multiply peaked 
clouds as the most evolved clouds. Most of the starless sources 
were in diffuse clouds, which are dark in the 24 jum MIPS im- 
ages. The prestellar sources were generally found in the vicin- 
ity of infrared objects, most likely more evolved YSOs. These 
YSOs are possibly responsible for introducing shocks into their 
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surrounding medium through stellar winds and outflows, thus 
supporting the formation of higher densities regions where grav- 
itational collapse can occur This implies that stars form close 
to other stars, which is a well-known fact and agrees with the 
observed clustering of newly formed stars (see, for example, 
lLada&Ladall200i . The limited number of both starless and 
prestellar sources in a total of 12 clumps, compared to the 
31 sources with a YSO indication, agrees with the findings of 
iMotte et aTl (|2007|) that the starless and prestellar phases are 
short and dynamic. 

3.4. CO depletion 

As a clump contracts, the level of de pletion of C-bearin; 
molecules like CO or CS increases (iKramer et alJ \199' 
[Bergin & Langer 1997J. The CO depletion in a clump can 
therefore be used as a time marker. O bservations of massive 
protostellar objects (P urcell et al.ll2009l) and clumps in IRDCs 
jPillai et al] 12007) show depletion of CO, indicating that these 
objects are akeady in a state of collapse and not in the chem- 
ically earliest stage. CO depletion as a time marker should be 
treated with some caution, because heating of the gas (by a YSO) 
can cause the CO to be desorbed to the gas phase, creating low 
CO-depletion levels as in the chemically earliest stages around 
the onset of collapse. 

We used J8o, a less common isotopologue of CO, to mea- 
sure the depletion in the clumps of the high extinction clouds. 
In contrast to CO, the line profile of J8o lines consisted of 
one Gaussian without any wings or self-absorption. This is an 
indication that the J8o(2-l) emission has a low or moderate 
optical depth, which would agree with previous observations 
that derived the J So optical depths by compa ring the ISo data 
with emission fr om the even rarer C'^O (iCrapsi et al.l l2005t 
iPillai et ani2007h . Therefore, in this work, we will assume that 
the J8o(2-l) line is optically thin. The J8o(2-l) line profiles 
were fitted with a Gaussian. Table IB. 81 lists our obtained inte- 
grated intensities, the LSR velocities, and the line widths. 

The column density of J8o was calculated in the optically 
thin limit (see Eq lA.4l i under the assumption of a beam fill- 
ing factor of unity, which is likely for a beam of 11" and 
clumps that have an average size of 20" in the 1.2 mm contin- 
uum. Second, we assume LTE and equate the excitation tem- 
perature to the rotational temperature, r,ot, derived from the 
NH3(1,1) and NH3(2,2) transitions. In very dense and cold re- 
gions, where NH3 inversion transitions are dominated by col- 
Hsions rather than by radiative excitations and de-excitations, 
Trot is a good measure for the k i netic temperature o f the system 
dWalmslev & Ung erechtslfTosl iDanbv et al.| [T988). The back- 
ground temperature was taken at 2.73 K. With the total column 
density of J 80, Nc'^o, known, we calculated the abundance of 
J80 relative to H2, Xo^o^ by comparing it to the hydrogen col- 
umn density: 

The canonical numbe r for;t'ciso in the nearby ISM is 1 .7 x 10"^ 
dFrerking et al.llT982l) . At larger distances towards the molecular 
ring and the Galactic center, this value is expected to be higher 
(Wilson & Rood 1994), which would increase the CO depletion 
we measured. The ratio between the observed abundance, Xobs, 
and the canonical value for the abundance, ;t'canonicai, then re- 
sulted in 77 the CO depletion factor: 




20 

,t(NH3) [K] 

Fig. 11. CO depletion versus the rotational temperature derived 
from NH3. The green crosses represent sources in diffuse clouds, 
the blue triangles clumps in peaked clouds, and the red squares 
sources in multiply peaked clouds. The black dashed line marks 
the trend of an increase in CO depletion with colder temperatures 
for all the clumps, while the green, blue, and red solid lines mark 
the trends for the sources in diffuse, peaked, and multiply peaked 
clouds. 



All derived quantities, A^c"*0' ^h,, Afc"*o> ^nd ^ ^e listed in 
Table |R9] In most of the clumps, J 80 appears to be depleted. 
We find a CO depletion of 77=0.5-5.5 with an average of 2.3, 
which is slig htly low er than what is found in clump s of IRDCs 
(77=2- 10, Filial et^ 2007). Obse rvations by, e.g., K ramer et al.l 
(119991) and Bacma nn et al T d2003h have shown that molecular de- 
pletion is enhanced at lower temperatures. In Fig.[TT] we plot for 
all clu mps the CO depletio n against the rotational temperature. 
While IKramer et alJ (1 19991) found an strong exponential depen- 
dence of CO depletion with temperature, our data are too scat- 
tered to differentiate between linear and exponential behavior 
Nevertheless we can see that higher 77 corresponds on average to 
colder sources (dashed black line). 

Figure[TT]shows that clumps in diffuse clouds (green crosses 
and green line) and peaked clouds (blue triangles and blue line) 
indeed follow the trend of higher CO depletion with colder tem- 
peratures. For clumps in multiply peaked clouds (red squares and 
red line), the situation is different; the clumps have a large spread 
in temperature as well as 77, and on average the opposite trend 
is seen: increasing 77 with higher temperature. This is mainly 
caused by a few clumps with high temperature and high 77, 
G017.19+0.81MM2, G014.63-00.63MM1, and G022.06+0.21 
MMl, all three of which have dust column densities above 

1023 



cm" and water masers. Possibly these sources were heated 
so fast by the embedded YSO that the CO was not freed into the 
gas phase yet and the depletion values are thus still high. 



4. Discussion 

4.1. Infall velocities 



In Sect. 13.11 we determined by binomial tests that the HCO^(l- 
0) line was the best suited of the molecular transitions that we 
had observed to study infall. For sources that had a double- 
peaked HCO^(l-O) line profile, indicating a systemic infalling 
motion, we calculated the infall velocity Vin according to 
iMvers et al.l(ll996l) as 



Vin = 



^'red ~ ^'blue 



/ 1 + e*^blue-7'dip)/7'dip 

ln( ^ „^ ), (9) 



A'canonical 
^obs 



(8) 



where Vred and Vbiue are the LSR velocities of the red and blue 
peaks, whose main beam temperatures are given by Tied and 
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Tbiue, respectively, and T^-p is the main beam temperature of the 
dip between the two peaks. All these quantities refer to the spec- 
tral shape of the optically thick line, except for Avthin- With the 
volume density of the clump rtu, and clump radius R, the in- 
fall velocity can be converted to a mass infall rate Min (after 
iMvers et al.|[l996h : 



Min = 4-nR nu^jim-iiVi, 



(10) 



where ju = 2.33 the mean molecular weight and mn the mass of 
an hydrogen atom. For the clump n^^ and clump radius R, we 
took the volume density and clump radius determined in Paper I 
from 1.2 mm dust continuum observations. 

For low-mass star formation, the mass infall rates onto a pro- 
tostar are between a few times ~ 10"^ Mq yr"' (iPalla & Stahlei 
1993[ iKirk et alj I200I and ~ IO ^'Mq yr ' ( Whit ney et al. 
19971) . This is too low to form high-mass stars, since their 
preste llar p hase is much shorter, less than 3 x 10"* yr, accord- 
ing to lMotte et all (|20()7|) . which would mean a mass infall rate 
of 3 X 10"^ Mq yr ' to form a 10, M© star. Theoretical mod- 
els predict a mass accretion onto a high- mass protostellar core 
(at radii < 0.1 pc) of Min of 10"^^ Mq yr"' dMcKee & Tanll2003l: 
[Banerjee & Pudritz 2007). Observations of massive clumps 
(with typical r adii > 0. 1 pc) firi d mass infall rates of 10 - 
10-4 Mg yr ' dFuller et al.1 12001 iPeretto et al.ll2006 l). In high- 
resolution observations of high-mass protostellar cores, the mass 
infall rate is indeed higher and covers the predicted order of mag- 
nitude (10^2- 10^4 Mq yr-', Beltran et al. 2006, 2008). Our mass 
infall results, given in Table |5j fall in the range of the observed 
values for massive star-forming clumps. 

To explore further the mass infall rate, we calculated the in- 
fall timescale using the clump mass obtained from the 1 .2 mm 
continuum (Paper 1): Tin = ^ yr. This timescale can be com- 
pared with the free-fall timescale - 3.3 x lO^n"'^^ yr, which 
was calculated using the density obtained from the same 1 .2 mm 
continuum observations (Paper 1). Both timescales are given in 
Table |5] In general, the free-fall timescale is longer than the 
estimated infall timescale, which indicates that the mass infall 
rate calculated in Eq. [10] is likely overestimated. The formula 
in Eq. [TO] assumes a spherically symmetric infall of the entire 
clump, while most likely the HCO^(l-O) emission is originat- 
ing in an optically thick low-density laye r in the outer parts of 
the clump (see also lLopez-Sepulcre et al.l (l2010 l)V We also cal- 
culated the HCO^(4-3) infall velocities and mass infall rates for 
the few sources that showed a double-peaked blue-excess line 
profile (Table |5]). Unfortunately, the error in the infall velocity 
of G024.94-00.15MM1 was quite high, leaving only G014.63- 
00.57 MMl to serve as a comparison. For the latter, the infall 
velocity derived from the higher excitation HCO^(4-3) line pro- 
file is half of the one derived from the HCO^(l-O) line pro- 
file, and the HCO+(4-3) infall timescale is twice the HCO+(1-0) 
timescale. At least for this one source, we find an indication that 
the HCO^(l-O) infall velocities and mass infall rates are indeed 
overestimated. 



4.2. Evolutionary sequence 

In Paper I, we defined three classes of high extinction clouds, 
based on the contrast of the clump to the cloud in the 1 .2 mm 
continuum and the number of clumps per cloud. These clouds 
comprise diffuse clouds, clouds with a low clump-to-cloud col- 
umn density contrast Cat^^ <2, peaked clouds, clouds with higher 
contrast 2<Cn^^ <3 having at least one clump, multiply peaked 
clouds, clouds with still higher contrast Ca;,,^ >3 that generally 



contained several clumps. Here, we review all the results in 
light of the hypothetical evolutionary sequence we proposed, 
i.e., from the quiescent diffuse clouds to the peaked ones, and 
finally to the multiply peaked clouds. 

4.2.1 . Diffuse clouds 

For the class of diffuse clouds, more than half of the investi- 
gated sources were not well-defined clumps but positions of dif- 
fuse emission within a cloud. The skewness parameter analysis 
found a slight excess of blue sources, indicating infall, but the 
binomial test showed that this result is not significant (see Table 
|6|l. The diffuse clouds are likely the phase where the clumps 
are beginning to collapse. This would agree with the clump 
masses found in Paper I, where on average we found less mas- 
sive clumps in the diffuse clouds than in the peaked or multi- 
ply peaked clouds (see also Table[T]l. Moreover, few YSO/hot 
core indications were found toward clumps in diffuse clouds, 
and half of the starless and prestellar sources were found toward 
the diffuse clouds (the other half was found toward the peaked 
clouds). A young evolutionary stage can also be identified by the 
CO depletion , since CO depletion requires a certain dens ity (> 
fewlO'* cm'l iKramer et al.lll999t iBergin & Tafallall2007h . Low 
depletion levels can indicate that the clump has not yet reached 
a high enough density. The mean depletion level of clumps in 
diffuse clouds (77=1.9) was slightly lower than that found toward 
the clumps in peaked and multiply p eaked clouds (?7=2.3). T he 
depletion timescale, ~ 5x 10^/»H,( Bergin & Tafal la"2007^. is 
~4 X lO"* yr for clumps in diffuse clouds, which is slightly longer 
than that for in peaked clouds (~3 x 10"^ yr) and multiply peaked 
clouds (~2 X 10'* yr). 



Table 6. Infall properties from the HCO^(l-O) line profile for 
different classes of clouds 



Class 




A'red 


^ skew 


A'.o, 


E 


P 


< 6v > 


Diffuse 


4 


2 


6 


8 


0.25 


0.34 


-0.001 


Peaked 


11 


2 


13 


17 


0.53 


0.01 


-0.63 


Multiply p. 


9 


8 


17 


20 


0.05 


0.50 


-0.72 



Notes. Columns are (from left to right) number of blue excess sources, 
number of red excess sources, total number of sources with significant 
HCO^(l-O) skewness, total number of sources, the excess parameter, 
the probability of the distribution to arise by chance, and the average 
HCO+(1-0) skewness. 



4.2.2. Peaked clouds 

In the peaked clouds, the distribution of blue and red sources 
from the HCO^(l-O) line profile significantly favors the blue 
sources with infall. The excess parameter was 0.53 with only a 
1% possibility for this distribution to arise by chance (Table |6]l. 
Apparently, most clumps in the peaked clouds are accreting ma- 
terial. The infall velocities and mass accretion rates for peaked 
(as well as the one for multiply peaked) clouds have a very wide 
range, whose mean lies far above the infall velocity/mass accre- 
tion rate of the diffuse clouds. All the SiO detections of clumps 
in peaked clouds coincide with sources that show infall based on 
the HCO^(l-O) skewness parameter (Fig.|9]l- The strong coiTe- 
lation between the infall characteristics and SiO outflows shows 
that accretion is accompanied by collimated outflows, which 
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Table 5. Infall velocities, mass infall rates, and timescales for the blue excess sources 



oource 




VinlO 


M . r. 


Tin 10 


l-m43 




1'in43 








(^Km s ) 


(,iu Mo yi ) 


UU yr) 


l^Km S ) 


(\ A/l .,,--1 \ 

(iu Mo yr ; 


(lu yi; 


yi; 


Diffuse clouds 


















G013.28-00.34.. 


. MMl 


1.09(0.15) 


1.5 


0.6 


- 


- 


- 


1.5 


G013.97-00.44.. 


. MMl 


0.31(0.13) 


- 


- 


- 


- 


- 


- 


Peaked clouds 


















G013.91-00.51.. 


. MMl 


0.69(0.10) 


2.2 


0.66 


- 


- 


- 


0.70 


G024.94-00.15.. 


. MMl 


1.34(0.12) 


3.3 


0.32 


0.34(0.66) 


0.85 


1.2 


0.74 




MM2 


3.00(0.13) 


6.2 


0.13 


- 


- 


- 


0.78 


G030.90+00.00. 


.. MM2 


1.76(0.15) 


9.4 


0.75 


- 


- 


- 


1.50 


G035. 49-00.30.. 


. MMl 


7.4(0.13) 


16.6 


0.06 


- 


- 


- 


0.78 




MM2 


84('0 1 0'l 


1.4 


0.5 








0.85 


G053.81-00.00.. 


. MMl 


7.4(0.10) 


9.6 


0.03 








0.55 


Multiply peaked 


clouds 
















G014.63-00.57... 


MMl 


2.7(0.14) 


22.0 


0.15 


1.30(0.35) 


10.6 


0.30 


0.38 




MM2 








0.90(0.28) 


2.7 


0.26 


0.38 


G017.19+00.81. 


.. MM2 


0.66(0.20) 


3.1 


0.46 








0.42 




MM3 


1.65(0.13) 


2.8 


0.26 








0.85 


G018.26-00.24.. 


MM3 


6.1(0.18) 


16.0 


0.08 








0.80 



Notes. Columns are (from left to right) source name, infall velocity; mass accretion rate; accretion timescale based on the HCO+(1-0) line, infall 
velocity; mass accretion rate; accretion timescale based on the HCO^(4—3) line, and free-fall timescale. 



leads to an enhancement of SiO. This correlation does not ex- 
ist for the earlier or later classes. 

4.2.3. Multiply peaked clouds 

Among the multiply peaked clouds, there were almost equal 
numbers of blue and red excess sources (based on the HCO^(l- 
0) skewness). These clouds seem to be hosting clumps in a vari- 
ety of evolutionary stages: from clumps with still ongoing accre- 
tion to clumps where the feedback from the formed young stars 
halts further accretion. 

Many clumps in these clouds presented SiO emission. Figure 
|9] shows that the number of clumps with SiO emission showing 
infall is about equal to those showing expansion (both infall and 
expansion characteristics are based on the HCO^(l-O) skewness 
parameter), meaning that the outflow continues even after large- 
scale infall has finished. 

The volume density of the multiply peaked clouds, nn, > 
10^ cm"^, is much larger than what is found for the peaked 
clouds. Multiply peaked clouds have the highest detection rate 
of the CH3CN(5^) K^Q and H2CO(4o3-3o4) lines, indicating 
that this class contains the most YSOs or hot molecular cores. 
In Paper I, we akeady found that the NH3 rotational temper- 
atures were on average higher for clumps in multiply peaked 
clouds than for the other classes (Table[T]). The K levels of the 
CH3CN(5^) line allowed a handful of rotational temperature 
estimates, ranging between 23 and 34 K. The increasing evi- 
dence for YSOs is likely the reason for the broadening of the 
line width found toward clumps in this class (observed for the 
N2H+(l-0), NH3(1,1), 18o(2-l), and Hi3cO+(l-0) lines) by in- 
creasing turbulence. 

Finally, three clumps in the multiply peaked clouds, 
G014.63-00.57 MMl, G017. 19+00.81 MM2, and 
G022. 06-1-00.21 MMl, all show various characteristics that 
are different to the properties of the bulk of the sample. 
Examples include their SiO intensity against column density, 
the CO depletion against temperature, and their N2H^(1- 
0)/N2H^(3-2) ratio. While these clumps have strong outflows, 
suggested by their SiO emission and the presence of water 



masers, they also have very high co lumn densities (theoreticall y 
enough to form massive stars, e.g.. lKruniholz & Mc Kee"2008h. 
They show strong CO depletion values while having elevated 
temperatures (~ 20 K) due to the rapid heating of the environ- 
ment by the YSO. Most likely, these three clumps are actively 
forming massive YSOs, which is why they have such extreme 
properties. 

5. Summary 

We performed spectral line observations of several molecules 
using the IRAM 30m and the APEX telescopes toward a sample 
of clumps in high extinction clouds spanning various evolution- 
ary stages. The skewness parameter was used to classify sources 
based on their line profiles as infall sources or outflow/expanding 
sources. The HCO^(l-O) and CO(3-2) lines were found to be 
more sensitive to infall and outflow than the HCO^(4-3) line. 
Only the HCO^(l-O) profiles showed a significant excess of in- 
fall sources, and therefore the HCO^(l-O) line profile was found 
to be most effective to investigate infall profiles. 

SiO emission, which is thought to be a tracer of shocks, was 
predominantly detected toward sources with infall. This finding 
is consistent with a scenario in which accretion is accompanied 
by collimated outflows that lead to an enhancement of SiO. 

Between the three cloud classes, the peaked clouds have 
a significant high excess of infall sources, while for the other 
classes no significant excess was found. Possibly, the infall phase 
has not yet begun in most diffuse clouds and already has come 
to an end in some of the clumps of the multiply peaked clouds. 
The infall velocities and mass infall rates suggest an increasing 
trend when going from the diffuse clouds (1.5 xlO"^ M© yr"') to 
the peaked (and multiply peaked) clouds (1. 5-16x10"^ yr"') 
suggesting that mass infall increases in these later stages. More 
measurements of clumps in diffuse clouds are, however, neces- 
sary to confirm this. 

From RADEX calculations, we found that the ratio of the 
N2H^(l-0) to N2H^(3-2) integrated intensity is sensitive to 
the hydrogen volume density, and together with the N2H+(1- 
0) optical depth can be used to derive good estimates of not 
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only the hydrogen volume density but also the intrinsic N2H^ 
column density. Apart from the ammonia rotational temper- 
atures, determined in Paper 1, we added several other tem- 
perature estimates. All estimates converge to a temperature 
range between 10 and 40 K, where the higher value is deter- 
mined by CH3CN, a molecule usually found in hot cores and 
UCHii regions. Similarly, various density estimates place the 
density in the order of 10^ cm"^. 

The presence of YSOs in the clumps of the high extinc- 
tion clouds was established by 24 //m emission, but also indi- 
rect evidence through outflows (SiO emission) and the pres- 
ence of hot molecular cores through CH3CN and H2CO emis- 
sion was gathered. We find that most of the YSOs are located in 
clumps of peaked and multiply peaked clouds, while almost no 
YSO indications were present in diffuse clouds. The presence of 
an YSO seems independent of the physical size of the clump, 
since YSOs are found over a large range of clump diameters. 
However, there seems to be (hydrogen) column density cutoff" at 
(4 - 5) X 10^^ cm"^, below which few or no YSOs are found. The 
cutoff is most clear in the 24 fim emission and in the presence of 
SiO emission. 

In summary, the diffuse clouds seem to contain clumps that 
are on the verge of infall. The first stages of star formation are 
ongoing in the peaked clouds, since a majority of the clumps 
is inf ailing and many clumps also show outffows. The multiply 
peaked clouds harbor the more evolved clumps, where infall has 
already stopped in some of them, and generally we find indi- 
cations of YSOs. Hence, the cloud morphology seems to have 
a direct connection with the evolutionary stage of the objects 
forming inside. 
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Appendix A: Formulas and equations used 

A.1. A/2H+ 

The excitation temperature Tex can be calculated as 

hv/k 



T - 

^ ex — 



(A.l) 



where A - Tmf{Jy{T) + Jv{Tbg) is the outcome of the CLASS hfs 
method fitting. Using ln{e -1- 1) - e, one can simplify the above 
equation as 

rex-Ar„/+l. (A.2) 

If we ignore the in the above equation, we find that Tex 
1 //. Since the average Tex was found to be around 10 K, ignoring 
the in the above equation would lead to an uncertainty of 
(on average) 10%. 

We calculated th e column density of N2H^, A^N2H+, following 
iBenson et al. 



AfN,H+ = 3.3 X 10" 



1 _ e-4.47/7-„ 



(A3) 
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where Av is the line width obtained from the hyperfine fit and 
Tex is the excitation temperature. 

A.2. CO 

In LTE, the total column density, Mot, of the CO molecule can 
be calculated from the integrated main beam brightness temper- 
ature of an optically thin line as follows: 

3h Q(rex) e^'/"^- /TmbAv 

Sn^S/ul gugKgnuc e""/*^^" - 1 /vC^ex) - /v(2.73) ' ^ ' ' 

where h Planck constant 

k the Boltzmann constant 

S Une strength 

Hd the dipole moment of the molecule 

Tex the excitation temperature 

QiTex) partition function at Tex, 

e(rex)^0.38(rex + 0.88) 

gu the rotational degeneracy for the upper energy 

gK the K degeneracy 

gnuc the nuclear degeneracy 

Eu the energy of the upper level 

V the transition frequency 

hv/k 



Jy{T) 



glty/kT_l 



The properties of the J8o J=2—>1 transition are: 5 = 
Jii/(2J,i + 1)=0.4, where /„ is the angular momentum quan- 
tum number of the upper J level, yL(£)=0. 1098 Debye=0. 1098 x 
10 '** e.s.u., gu=2Ju + 1=5, gK=l, gnuc=l, Eu=kx 15.81 erg, and 
v=219.56GHz. 

The general formula for the CO column density, valid in the 
optically thick and thin regime is 

Mo. = . , 7 tAv, (A.5) 



Sn^Sfll gugKgnuc e'"'/*^' - „ 
where all the quantities are as defined above. 
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Fig. 6. N2H^(l-0) optical depth versus the ratio of N2H^(l-0) and N2H^(3-2) integrated intensities. The four plots represent 
RADEX calculations with different temperatures: top left 10 K; top right 15 K; bottom left 20 K; bottom right 25 K. The black 
dashed contours represent the logarithm of the line width times the column density, log(NN2H+ / Av) with a contour step of 0. 1 . The 
red dashed contours are logina-,) with a contour step of 0.2. The lowest and highest contour values are given in each plot. The values 
of the sources in Table|4]are plotted by black dots. 
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Table B.l. Molecular line survey overview: ' V for a detection and when not observed. For non-detections 3cr upper limits are given in Kelvin. Positions are expressed in J2000 coordinates. Water maser 

detections from Paper I. 
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MMl 


19 


23:12.4 


15:13:35 


V 


V 




V 




<0.18 


<0.22 




V 


V 




n.d. 




MM2 


19 


23:09.2 


15:12:42 


V 


V 




V 




<0.20 


<0.26 




V 


V 






G053.81-00.00... 


MMl 


19 


30:55.7 


18:29:55 


V 


V 




V 


V 


V 


<0.22 


V 


V 


V 




n.d. 



Table B.2. Infall and outflow study using the HCO*(1-0), C0(2-l), HCO*(4-3) and CO(3-2) lines. When the source was not observed the field is marked as 



Source name 


HCO+(1-0) 


CO(2-l) 


HCO+(4-3) 


CO(3-2) 




1 reel 


Profile Profile Wings/Range 


Profile 6v 


Profile 6v 




/(km s-') 







G012.73-00.58... 


MMl 


- 


0.19(0.01) 


- 


- 


y/11 












MM2 


0.43(0.04) 


0.51(0.13) 


R 


- 


y/13 










G013.28-00.34... 


MMl 


1.58(0.34) 


-0.59(0.13) 


B 


B 


y/14 










G013.91-00.51... 


MMl 


1.63(0.17) 


-0.64(0.19) 


B 


R 


y/27 


B 


-0.46(0.21) 


- 


- 


G013.97-00.45... 


MMl 


1.24(0.17) 


-0.73(0.13) 


B 


- 


n 










G014.39-00.75A.. 


MMl 
MM2 


- 


-0.52(0.03) 
-0.45(0.03) 


B 
B 


B 
B 


y/18 
y/13 










G014.39-00.75B.. 


MM3 


- 


-5.79(0.29) 


B 


- 


y/15 










G014.63-00.57... 


MMl 


1.05(0.09) 


- 


- 


B 


y/30 


B 


-0.59(0.21) 


B 


-2.06(0.21) 




MM2 


0.45(0.13) 


1.18(0.13) 


R 


B 


y/23 


B 


-0.78(0.21) 


- 


- 




MM3 


- 


0.55(0.04) 


R 


B 


y/16 












MM4 


- 


0.22(0.03) 


- 


- 


y/12 










G016.93+00.24... 


MMl 


- 


-0.62(0.13) 


B 


- 


n 










G017.19+00.81... 


MMl 


0.35(0.04) 


0.32(0.13) 


R 


R 


y/29 


R 


0.37(0.21) 


R 


1.07(0.21) 




MM2 


1.47(0.15) 


-0.48(0.13) 


B 


R 


y/43 


R 


0.55(0.21) 


R 


1.15(0.21) 




MM3 


2.39(0.36) 


-0.35(0.13) 


B 


B 


y/35 












MM4 


0.40(0.08) 


0.27(0.16) 


R 


R 


y/14 


- 


-0.06(0.23) 


R 


0.48(0.23) 


G018.26-00.24... 


MMl 


- 


-1.18(0.04) 


B 


B 


y/20 


R 


0.42(0.22) 


B 


-1.45(0.22) 




MM2 


- 


-1.03(0.14) 


B 


B 


y/12 


- 


-0.11(0.22) 


B 


-0.98(0.22) 




MM3 


1.68(0.65) 


-1.12(0.14) 


B 


B 


y/25 


- 


0.00(0.22) 


B 


-1.18(0.22) 




MM4 


0.89(0.35) 


- 


- 


R 


y/23 


- 


-0.04(0.23) 


B 


-1.14(0.23) 




MM5 


- 


0.49(0.16) 


R 


B 


y/15 










G022.06+00.21... 


MMl 


0.23(0.07) 


0.79(0.14) 


R 


R 


y/46 


- 


-0.11(0.22) 


R 


1.31(0.22) 




MM2 


0.55(0.21) 


0.76(0.15) 


R 


R 


y/12 










G023.28-00.12... 


MMl 


- 


-0.47(0.14) 


B 


- 


n 


- 


-0.23(0.22) 


R 


0.87(0.22) 




MM2 


0.26(0.07) 


0.33(0.22) 


R 


- 


y/14 










G024.37-00.15... 


MMl 




-0.59(0.22) 


B 


B 


y/11 












MM2 




-1.11(0.26) 


B 


R 


y/17 


- 


-0.15(0.31) 


B 


-1.11(0.31) 


G024.61-00.33... 


MMl 


- 


-0.72(0.13) 


B 


B 


y/27 


R 


0.26(0.22) 


B 


-1.75(0.21) 




MM2 


- 


-0.07(0.14) 


- 


B 


y/40 










G024.94-00.15... 


MMl 


2.03(0.24) 


-0.80(0.14) 


B 


B 


y/21 


B 


-0.43(0.28) 


B 


-0.99(0.22) 




MM2 


2.24(0.30) 


-0.66(0.16) 


B 


B 


y/27 










G025.79+00.81... 


MMl 


0.44(0.10) 


1.08(0.14) 


R 




y/10 












MM2 




0.25(0.06) 


R 




n 










G030.90+00.00A.. 


MMl 




0.66(0.03) 


R 


R 


y/20 




-0.04(0.16) 


R 


0.63(0.22) 


G030.90+00.00B.. 


MM2 


1.90(0.24) 


-0.42(0.18) 


B 


B 


y/13 




-0.17(0.26) 




-0.19(0.24) 


G030.90+00.00C.. 


MM3 




-0.27(0.02) 


B 


B 


y/26 










G034.71-00.63... 


MMl 




-0.11(0.01) 






y/43 












MM2 




-0.35(0.07) 


B 


R 


y/25 


B 


-0.40(0.22) 


R 


0.35(0.22) 




MM3 




1.09(0.04) 


R 


R 


y/22 










G034.77-00.81... 


MMl 




-0.76(0.07) 


B 


B 


n 











Table B.2. continued. 



kjtJUlcC llalllC 


















— Jw — 


6v 


Profile 


Profile 


Wings/Range 
/(km s-i) 


Profile 6v 


Profile 6v 


nO^zL 8 S -1-00 zL^ 


Mivn 

iVliVl i 






I\ 


D 


n 






nO^S 4Q-00 ^OA 


iVliVl 1 






1 ) 


Ti 

1 ) 


v/1 f\ 

y/io 








MM? 








13 
1 ) 


v/1 1 








MM3 




-0.18(0.01) 


- 


B 


y/U 






G037.44+00.14A.. 


MMl 










n 






G037.44+00.14B.. 


MM2 




0.00(0.20) 






n 






G050.06+00.06... 


MMl 


0.97(0.12) 


0.0(0.14) 






y/14 








MM2 




0.00(0.01) 






y/17 






G053.81-00.00... 


MMl 


3.02(0.33) 


-0.34(0.13) 


B 


B 


y/34 







Notes. Columns are (from left to right) source name, ratio of the blue over red peak intensity of the HCO'^(l-O) line, HCO*(1-0) skewness parameter, HCO*(1-0) profile, CO(2-l) profile, CO(2-l 
wings/CO range of wings, HCO*(4-3) profile, HCO*(4-3) skewness parameter, CO(3-2) profile, and CO(3-2) skewness parameter. 



Table B.3. Molecular line survey results for HCO*(1-0) and H"CO*(1-0); columns 7-1 1 are based on Gaussian fits. 



Source name 








HCO+(1-0) 








H'^CO+(1-0) 








T'blue 


VLSR,blue 


^red 


Vl CD 


r 


vlsr 


r 


Vlsr 


AV 






(K) 


(kms-1) 




(km s"') 


(K) 


(km s-i) 


(K) 


(kms-i) 


(km s"' 


/^mo T2 r\r\ 
LrUlZ. / J— UU.JO... 


MMl 






- 


- 


1 TO/'n (\ci\ 

1. /z(u.uy) 


o.oz(U.Ul; 






A ^Qtr\ A^\ 




MMZ 


l.U4(U.UDj 


3.iD(U.Zj) 


2.41(0.06) 


6.87(0.25) 






i.iy(u. 11) 


^ CA/A A0\ 


A TO/A AC^ 

U. /Z{().()j ) 


/~'r\i^ oc r\r\ q/i 
CjU 1 J .ZS— UU. j4. . . 


TV /TA/ri 

MM i 


l.zU(U.iU) 


/I A 1 T/A 
4U.i /(U.zj) 


0.76(0.10) 


42.49(0.25) 






y) X)v>) 


A 1 '2A/'A AA\ 


1 O 1 /A 1 <\ 

i.yi(u.ij) 


/^rvi Q o 1 r\r\ ^ 1 


A TA /T 1 
MM 1 


Z.Uo(,U.Uo; 


ZZ.Uo(U.Zj ) 


1.27(0.08) 


24.05(0.25) 






1 n{C\ c\n\ 
1. /(U.U/) 


OO m/A AO\ 

zz.y /(U.UZ) 


1 QO/'A r\A\ 

i.jy(u.u4) 


/^niQ m r\r\ 
CjUl J.y /— UU.4D... 


A yfA yf 1 

MMl 


1 "^/T* (\(\\ 

1.4j(U.Uy; 


1 Q < 1 /"A 


1.15(0.09) 


20.88(0.25) 






U.oz(U.Uo) 


OA A'2/'A A^\ 


O AO /A 1 1 \ 

Z.Uo(U.l 1) 


LrU14.jy— UU. /jA.. 


A >fA >f 1 

MMl 






- 


- 


z.jy(u.iz) 


1 /.ZU(U.U3 ) 


1 no/n r\o\ 
l.UZ(U.Uo ) 




1 OO/A A7\ 

l.ZZ(U.U/) 




A JA JO 

MMz 






- 


- 


U.oo(U.lz) 


lo.yy(U.Uo) 


A TA/A A/;\ 


1 T CO/A A'3\ 

1 /.Jo(U.U3) 


1 O/A AO\ 

l.J(U.Uo) 


(jU14.oj— UU.J /... 


A JA J1 

MMl 


j.io(U.l j) 


iD.iy(U.Zj) 


3.04(0.13) 


20.30(0.25) 






1.Zj(U.Uj ) 


lo.jj(U.U3 ) 


A </A A7\ 

Z.43(U.U/ j 




MMZ 


U.jo(,U.l 1 ; 


17 1 A^(\ 

1 / .14(U.Zj ) 


1.24(0.11) 


20.28(0.25) 






1.j4(U.Uj ) 


-IQAA /'A AO\ 

lo.44(U.Uz; 


1 <A/A A/l\ 

l.JO(l).U4) 




A /fA yfo 

MM J 






- 


- 


l.ZU(U.Uy) 


lo.z4(U.Uo) 


U.Ol(U.Uo) 


in Al /"A 

1 / .41(U.Uj) 


1 ^A/A 1 0\ 




MM4 






- 


- 


U.jj(U.U/) 


iy.4z(u.uy) 




iy.z4((J.(J4) 


A OA/A 1 0\ 

U.oU(U.lZ) 


(jU10.yj+UU.Z4... 


MM 1 


z.jj(U. i i) 


T2 1 /I /A 0<\ 

z3.i4(U.Zj) 


- 


- 






1.1 3(U.Uj) 


Ol "7"7/A A1 \ 
z3. / /(U.Ul ) 


1 A1 /A A/1\ 

i.Ui(U.U4) 


LrUl /.ly+UU.ol... 


\ /I \ /I 1 

MM 1 


i. /0(U. 14) 


zi.&o(U.Zj ) 


5.05(0.14) 


25.48(0.25) 






l.Z3(U.Uo) 


0< AA/'A AO\ 


1 <A/A A<\ 




A /TA yfO 
MMZ 


4. / /(U.ZU) 


zi.yj^u.zj ) 


3.25(0.20) 


24.82(0.25) 






1.jj(,U.Uj J 


OO '7Q/'A AO\ 

zz. /yi^u.uzj 


1 70/A A^^ 




A /fA yfo 

MM J 


4.4/(^1). Iz) 


ZZ.l6(U.ZD) 


1.86(0.12) 


24.44(0.25) 






I. /o(U.Uo) 


OO ^'2^n ao\ 
ZZ./j(U.Uz) 


l.j4(U.U4) 




MM4 


U.yj^U.Uj; 


ly.oy(u.zD) 


2.33(0.05) 


22.69(0.25) 








OO 1 /'A A^\ 

ZZ. 1(1). Uj ) 


O 1 O/A 1 1 \ 

z. iy(u.i 1 ) 


vjU 1 o .Zo— UU . Z4 . . . 


TV /TIV A 1 

MM 1 






- 


- 


1.j1(U.1 j 




A "7 O/A A/1 \ 

U. /o(U.U4) 


AC /I /A A/1 \ 


O /I O/A AO\ 

z.4y(u.uy) 




MMZ 


l.JD(U.iZj 


oj.jy(u.zj ) 


- 


- 






1 AA/A A<\ 


/^T TO/'A A'J\ 

/. /y(u.u3 ) 


O 1 /I /A A0\ 

Z. 14(U.Uo ) 




MM J 




/:c A/I /'A 
0J.04(U.ZD) 


0.69(0.17) 


71.45(0.25) 






A QO/'A A/l\ 


AO/'A A/l\ 


O A7/A AO\ 

z.o/(u.uy) 




MM4 


r\ QAtf\ 1 n\ 

U.64(U.l /; 


OJ.ZO(U.ZD) 


0.94(0.17) 


71.66(0.25) 






A TO/'A A/l\ 


zTQ Q-2/A A<\ 


o n</A 1 1 \ 
Z.yD(U.l 1) 




MMj 


ATI /A O C\ 
U. / l(U.Z5) 


0/.Za(U.z5) 


- 


- 






A TC/A A/1\ 

U./j(U.U4) 


0D.J0(U.U4) 


O /t 1 /A 1 1 \ 

Z.43(U.l 1) 


OUZZ.Uo+UU.Zi... 


MM 1 


A C/i/'A OA\ 


4o.j4(U.Zj ) 


3.81(0.20) 


52.95(0.25) 






1 '20/'A A/^\ 

1.jZ(U.U0J 


C 1 Q C/A A'2\ 


O A1 /A A0\ 

Z.Ul(U.Uo) 




MMZ 


A ylO/'A 1 A\ 

U.4z(U.lU; 


4y.jD(U.ZD) 


0.77(0.10) 


53.16(0.25) 






A /^^/A A^\ 


C 1 OO/'A A<\ 


1 7</A 1 \ 
1. /j(U.l ) 


LrUZj.Zo— UU.IZ... 


A yTA/TI 

MMl 






- 


- 


j.zj(U.lz) 




A nn{f\ c\z\ 
U. / /(U.Uj j 


on Al/A A/l \ 

yy.U3(u.U4) 


O AQ/A AO\ 




A yfA 

MMZ 


A /lO/'A AO\ 

u.4z(u.uy; 


yo.i j(u.zj ) 


1.63(0.09) 


99.94(0.25) 






A QO/'A AA\ 


OO ^O/'A AO\ 

yy.jz(u.uo) 


1 '2 1 /A 1 0\ 

l.jl(U.lo) 


(jUz4.j /— UU.Ij... 


A JA J1 

MMl 


ATI /"A A/^\ 


J /. / /(U.ZD) 


- 


- 






A 'JO/'A A/1 \ 


CO 0/1 /'A A0\ 

jy.z4(u.Uo ) 


O /I O/A 1 7\ 

z.4y(u.i / ) 




A /TA TO 

MMZ 






- 


- 


f\ '2'7/A (\Q\ 


CO T/A 


A /I ^/A A/^\ 


OQ/A A0\ 


O 07/A 00\ 


(jUz4.01— UU.JJ... 


A TA yfl 
MMl 


l.Uo(U.14j 


41.5o(U.Zj) 


- 


- 






1 /I /I /'A A/^\ 

1.44(U.Uo; 


/lO '70/'A AO\ 

4z. /y(u.uz) 


1 70/A A^\ 

1. /Z(U.Uj ) 




MMZ 


1 OQ/A AO\ 


4j.jZ(U.Zj ) 


- 


- 






A TA/'A AA\ 


A Q. /^I/'A A/1 \ 


1/11 /A AO\ 

i.4i(u.uy) 


UUZ4 .y 4-UU .15... 


A JA 

MMl 


O iCO/A 1 A\ 
Z.OZ(U.iU) 


45. / /(U.Z5) 


1.29(0.10) 


49.28(0.25) 






A OO/A AC\ 

U.oo(U.Uj) 


/I T OT/A A/1 \ 

4/.3 /(U.U4) 


O AA/A AO\ 

Z.UU(U.Uo) 




A /TA /TO 

MMz 


o A*::/A An\ 

z.uo(u.uy) 


4o.y3(U.Z5) 


0.93(0.09) 


49.91(0.25) 






A C /I /A A/1 \ 

U.j4(U.U4) 


/I O 1 T/A AC\ 

45.1 /(U.Uj) 


1 07/A 1 1 \ 
1.0 /(U.l 1) 


(jUzd. /y+uu.oi... 


A/TA/T 1 
MMl 


A /lO/'A A0\ 


4y.l6(U.Zj ) 


1.11(0.08) 


51.57(0.25) 






A AO/'A A^\ 


/lO Q^/'A A/l\ 

4y.6j(U.U4) 


1 AA/A AO\ 

i.ou(u.uy) 




AJA /fO 

MMZ 






- 


- 


U.zj(U.U4) 


jU.z/(U.14) 


A O/i/'A AQ\ 

U.zo(U.Uo) 


CA n'2^n 1 1 \ 


u.y /(u.z4) 


CjU jU.yU+UU.UUA. . 


A TA J1 

MMl 






- 


- 


Z.Ol(U.14) 


/o.oo(U.Uj) 


A OO/'A AA\ 


O^/'A A/I\ 


O 1 /I /A AO\ 

z.i4(u.uy) 


CjUiU.yU+UU.UUr}.. 


MMZ 


/I I/A 1 1 \ 

z.43(U.l 1 ) 


1 <0/A A<\ 

yi.jy(u.ujj 


1.28(0.11) 


95.06(0.25) 






A TO/A A<\ 


OO 1 /'A A/^\ 

yz.oi(u.uo) 


O 0/1 /A 1 '2\ 

z.y4(u.ij) 


cjUju.yu+uu.uuc. 


MM 3 










yj. /4(U.U3 ) 


y3.yi(u.iz) 


A <0/'A A^\ 


0/1 'J'J/'A A/1 \ 

y4.33(U.U4) 


1 <0/A 1 0\ 
l.JO(U.iZ) 


oUj4. / 1— UU.Oj... 


IVllVl 1 










Z.joyyj.zj ) 


44.jZ(^U.U / J 


1.UU(,U.Uj J 


44. /o^U.Uj J 


OA/A AQ^ 
Z.ZOt^U.Uo) 




MM2 


2.11(0.11) 


45.23(0.25) 










0.76(0.04) 


45.83(0.03) 


1.73(0.08) 




MM3 


0.64(0.07) 


48.41(0.25) 










0.88(0.05) 


46.49(0.03) 


1.76(0.07) 


G034.77-00.81... 


MMl 










0.59(0.07) 


42(0.07) 


0.23(0.03) 


43.68(0.1) 


2.21(0.21) 


G034.85+00.43... 


MMl 










0.83(0.07) 


56(0.04) 


0.53(0.06) 


55.72(0.03) 


0.90(0.08) 


G035.49-00.30A.. 


MMl 


1.14(0.05) 


54.06(0.25) 


0.41(0.05) 


57.7(0.25) 






0.46(0.04) 


55.61(0.07) 


2.17(0.15) 



Table B.3. continued. 



Source name 








HCO-'(l-O) 








X T 1 ^ 4-/1 \ 

H'^CO^(l-O) 








T'blue 


VLSR,blue 


Tred 


VLSR.red 


Tmb 


vlsr 


Tmb 


Vlsr 


AV 






(K) 


(km s-i) 


(K) 


(km s-1) 


(K) 


(km s-i) 


(K) 


(kms-i) 


(km s"' 




MM? 




4.4. MfTl 7S~> 


1 ofii'n n6~> 


46 7A(0 75~> 








4S ^Qf'O Ci^\ 
,jy\\j,\jj ) 


1 f\H(C\ C\l\ 

1 .LtOV^L/.L/ / ) 




MM3 


2.92(0.07) 


45.02(0.25) 






_ 


_ 


0.92(0.05) 


45.85(0.04) 


1.99(0.09) 


G037.44+00.14A.. 


MMl 














0.93(0.07) 


18.26(0.02) 


0.64(0.04) 


G037.44+00.14B.. 


MM2 










0.27(0.17) 


40.46(0.08) 


0.28(0.17) 


40.46(0.08) 


1.33(0.17) 


G050.06+00.06... 


MMl 


0.98(0.06) 


53.23(0.25) 


1.01(0.06) 


55.36(0.25) 






0.62(0.04) 


54.44(0.04) 


1.89(0.09) 




MM2 










2.96(0.06) 


54.93(0.01) 


0.60(0.05) 


54.93(0.05) 


1.73(0.11) 


G053.81-00.00... 


MMl 


3.66(0.10) 


23.69(0.25) 


1.21(0.10) 


25.74(0.25) 






0.76(0.04) 


24.13(0.02) 


1.29(0.06) 



Notes. Columns are (from left to right) HCO*(1-0) blue peak intensity and its LSR velocity, HCO"^(1-0) red peak intensity and its LSR velocity, HCO*(1-0) peak intensity and its LRS velocity 
H"CO^(1-0) peak intensity and it LSR velocity, FWHP of the H"CO-"(1-0) line. » 
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Table B.4. Properties of N2H'^(l-0) based on Gaussian fits. 



Source name 




J TmbAv 


Vlsr 


Av 


Ttot 


Tex 


7Vn2H+ 


X 






(Kkms"') 


(km s""') 


(km s"') 




(K) 


(10^2 cm-2) 


(10-') 


G012.73-00.58... 


MMl 


3.7(0.1) 


6.68(0.01) 


0.42(0.01) 


4.11(1.18) 


5.3(1.6) 


5.3(2.2) 


0.26 




MM2 


6.7(0.2) 


6.35(0.01) 


0.83(0.04) 


5.97(0.91) 


4.6(0.8) 


12.1(2.8) 


- 


G013. 28-00.34.. 


MMl 


9.3(0.4) 


41.33(0.02) 


1.68(0.07) 


5.58(0.84) 


4.3(0.8) 


20.6(5.0) 


4.0 


G013.91-00.51... 


MMl 


12.9(0.4) 


22.75(0.01) 


1.19(0.03) 


5.34(0.45) 


5.7(0.6) 


22.0(3.0) 


1.0 


G013.97-00.45... 


MMl 


8.5(0.3) 


19.87(0.03) 


2.07(0.09) 


1.25(0.47) 


7.8(3.0) 


15.2(8.2) 


- 


G014.39-00.75A.. 


MMl 


6.6(0.3) 


17.64(0.02) 


1.14(0.07) 


0.78(0.76) 


- 


- 


- 




MM2 


6.9(0.3) 


17.12(0.02) 


0.92(0.05) 


0.82(0.86) 


- 


- 


- 


G014.39-00.75B.. 


MM3 


6.7(0.3) 


21.15(0.01) 


0.72(0.04) 


8.99(1.55) 


4.4(0.9) 


14.8(4.0) 


1.0 


G014.63-00.57... 


MMl 


45.8(2.3) 


18.45(0.01) 


2.46(0.02) 


2.35(0.13) 


13.3(0.8) 


89.0(7.4) 


0.52 




MM2 


18.0(0.7) 


18.23(0.01) 


1.27(0.03) 


3.90(0.42) 


7.6(0.9) 


27.8(4.5) 


0.5 




MM3 


7.7(0.3) 


17.44(0.03) 


1.34(0.08) 


3.71(0.93) 


4.8(1.3) 


13.0(4.8) 


1.4 




MM4 


2.8(0.2) 


19.03(0.03) 


0.63(0.06) 


6.66(2.54) 


3.7(1.6) 


7.3(4.2) 


2.4 


GO 16.93 +00.24... 


MMl 


8.7(0.3) 


23.63(0.02) 


0.84(0.05) 


4.17(1.05) 


5.9(1.6) 


12.8(4.8) 


1.0 


G017.19+00.81... 


MMl 


21.0(0.6) 


24.83(0.01) 


1.54(0.04) 


0.43(0.37) 


- 


- 


- 




MM2 


32.9(1.1) 


22.65(0.01) 


2.21(0.04) 


1.0(0.23) 


27.7(6.4) 


135.0(43.8) 


- 




MM3 


25.1(0.7) 


22.50(0.01) 


1.74(0.02) 


0.94(0.17) 


29.(0.3) 


109.7(28.0) 


- 




MM4 


4.5(0.2) 


22.00(0.03) 


2.10(0.05) 


0.44(0.45) 


- 


- 


- 


G018.26-00.24... 


MMl 


23.5(1.0) 


68.45(0.02) 


2.50(0.06) 


2.10(0.32) 


8.7(1.4) 


37.4(8.3) 


1.1 




MM2 


26.0(1.0) 


67.87(0.02) 


2.16(0.01) 


4.00(0.22) 


7.1(0.4) 


43.2(3.4) 


1.7 




MM3 


37.4(1.3) 


68.65(0.01) 


2.70(0.01) 


0.18(0.01) 


- 


- 


- 




MM4 


30.6(1.0) 


68.72(0.01) 


2.94(0.02) 


0.18(0.21) 


- 


- 


- 




MM5 


16.8(0.8) 


66.24(0.02) 


1.75(0.06) 


3.47(0.45) 


6.4(0.9) 


25.5(5.0) 


1.6 


G022.06+00.21... 


MMl 


28.7(1.3) 


51.25(0.01) 


2.50(0.01) 


1.17(0.01) 


17.7(0.1) 


76.6(0.5) 


0.48 




MM2 


8.3(0.4) 


51.59(0.02) 


1.74(0.08) 


1.57(0.52) 


7(2.4) 


13.4(6.4) 


0.52 


G023.28-00.12... 


MMl 


17.9(0.6) 


98.87(0.02) 


2.11(0.06) 


0.45(0.32) 


- 


- 


- 




MM2 


4.07(0.2) 


99.35(0.03) 


0.95(0.09) 


4.95(1.98) 


3.9(1.8) 


8.9(5.5) 


- 


G024.37-00.15... 


MMl 


13.5(0.8) 


58.87(0.02) 


2.53(0.05) 


1.33(0.27) 


8.6(1.8) 


23.5(6.9) 


1.4 




MM2 


12.8(0.5) 


56.18(0.02) 


2.28(0.07) 


1.55(0.32) 


7.9(1.7) 


21.4(6.4) 


1.3 


G024.61-00.33... 


MMl 


16.2(0.5) 


42.82(0.01) 


2.12(0.04) 


0.19(0.54) 


- 


- 


- 




MM2 


10.6(0.4) 


43.56(0.01) 


1.66(0.04) 


0.18(0.06) 


- 


- 


- 


G024.94-00.15... 


MMl 


20.0(0.8) 


47.24(0.02) 


2.48(0.05) 


1.30(0.20) 


12.1(1.8) 


41.6(8.9) 


0.88 




MM2 


11.9(0.4) 


48.22(0.01) 


1.92(0.05) 


0.36(0.30) 


- 


- 


- 


G025.79+00.81... 


MMl 


11.1(0.3) 


49.78(0.01) 


1.56(0.05) 


1.35(0.34) 


10.4(2.9) 


20.7(8.2) 


- 




MM2 


3.0(0.2) 


49.78(0.06) 


1.90(0.24) 


1.62(0.91) 


- 


- 


- 


G030.90+00.00A.. 


MMl 


22.8(0.8) 


75.21(0.01) 


2.38(0.02) 


0.18(0.10) 


- 


- 


- 


G030.90+00.00B.. 


MM2 


25.2(0.3) 


92.95(0.01) 


3.15(0.02) 


0.26(0.12) 


- 


- 


- 


G030.90+00.00C.. 


MM3 


13.6(0.6) 


94.29(0.02) 


1.50(0.06) 


1.76(0.51) 


9.6(2.9) 


22.4(9.4) 


1.1 


G034.71-00.63... 


MMl 


11.7(0.7) 


44.57(0.04) 


2.31(0.10) 


2.61(0.61) 


5.3(1.4) 


18.5(6.6) 


1.7 




MM2 


11.5(0.5) 


45.54(0.02) 


1.58(0.08) 


4.90(0.91) 


4.9(1.1) 


20.9(6.2) 


1.5 




MM3 


8.9(0.3) 


46.24(0.01) 


1.59(0.05) 


1.47(0.35) 


8.1(2) 


14.7(5.1) 


1.7 


G034.77-00.81... 


MMl 


1.6(0.1) 


43.29(0.06) 


1.60(0.19) 


1.11(1.26) 


- 


- 


- 


G034.85+00.43... 


MMl 


3.5(0.1) 


55.68(0.02) 


0.72(0.05) 


8.20(2.04) 


3.6(1.1) 


9.8(3.9) 


- 


G035.49-00.30A.. 


MMl 


8.3(0.3) 


55.28(0.01) 


2.14(0.04) 


0.18(0.02) 








G035.49-00.30B.. 


MM2 


9.4(0.4) 


45.40(0.01) 


1.37(0.04) 


5.52(0.64) 


4.6(0.6) 


18.4(3.3) 


1.3 




MM3 


14.8(0.5) 


45.69(0.01) 


1.86(0.04) 


3.28(0.30) 


6(0.6) 


23.0(3.2) 


1.7 


G037.44+00.14A.. 


MMl 


3.4(0.1) 


18.14(0.01) 


0.57(0.02) 


4.19(0.90) 


4.5(1) 


5.6(1.7) 




G037.44+00.14B.. 


MM2 


1.0(0.1) 


40.15(0.13) 


1.07(0.35) 


0.92(5.20) 








G050.06+00.06... 


MMl 


5.3(0.2) 


54.23(0.02) 


1.32(0.07) 


3.83(0.82) 


4.1(1) 


10.3(3.4) 


1.4 




MM2 


4.3(0.2) 


54.93(0.03) 


1.32(0.08) 


1.85(0.85) 


5(2.4) 


6.8(4.5) 


0.63 


G053.81-00.00... 


MMl 


9.0(0.3) 


24.11(0.01) 


1.64(0.03) 


0.45(0.25) 









Notes. Columns are (from left to right) source name, n2hp(l-0) peak main beam brightness temperature, LSR velocity at the N2H*(l-0) peak 
intensity, FWHP HCO+(1-0) line width, N2H+(l-0) total optical depth, N2H+(l-0) excitation temperature, N2H+ column density, and the N2H+ 
abundance relative to the H2 column density obtained from the 1.2 mm observations (Paper I). 



Table B.5. Observed properties of N2H*(3-2) and HCO*(4-3). When the source was not observed the field is marked as ' 



Source name 






N2H+(3-2) 








HCO+(4-3) 














1 7\jTidv 


"lsr 


Av 


J -* MB*^'^ 


"lsr 


Av 


^ blue 




red 


*^LSR,red 








(K Vm 




^^Mll a ) 


(K Vm s:"^^ 


(\crr\ 


\Kni S ) 












G013.91-00.51.. 


MMl 


2.4(0.2) 


22.77(0.08) 


1.56(0.2) 


4.2(0.4) 


22.43(0.06) 


1.38(0.16) 


- 


- 


- 


- 




G014.39-00.75A. 


MMl 


1.9(0.3) 


17.67(0.1) 


1.32(0.26) 


















G014.63-00.57.. 


MMl 


25.6(0.4) 


18.63(0.04) 


5.34(0.08) 


27.9(0.4) 


17.75(0.04) 


5.05(0.1) 


6.92(0.30) 


17.1(0.2) 


4.68(0.30) 


20.4(0.2) 






MM2 


8.2(0.4) 


18.42(0.07) 


3.31(0.2) 


7.6(0.3) 


18.29(0.08) 


3.75(0.19) 


1.99(0.26) 


17.3(0.2) 


1.78(0.26) 


18.9(0.2) 




G017.19+00.81.. 


MMl 


9.8(0.4) 


25.18(0.05) 


2.88(0.14) 


8.8(0.5) 


25.05(0.06) 


2.31(0.14) 


- 


- 


3.50(0.45) 


25.5(0.2) 






MM2 


18.7(0.4) 


22.82(0.04) 


4.21(0.1) 


25.5(0.7) 


22.78(0.08) 


5.99(0.19) 


3.92(0.45) 


21.4(0.2) 


5.26(0.45) 


23.8(0.2) 


,^ 




MM4 


1.7(0.3) 


21.57(0.15) 


1.88(0.37) 


4.1(0.5) 


21.97(0.15) 


2.61(0.36) 


_ 


- 


1.4(0.44) 


22.7(0.2) 


r 


G018.26-00.24.. 


MMl 


8.6(0.3) 


68.52(0.10) 


5.05(0.21) 


9.2(0.7) 


68.67(0.21) 


5.38(0.48) 


1.55(0.46) 


66.9(0.2) 


2.12(0.46) 


69.4(0.2) 






MM2 


7.6(0.4) 


67.95(0.10) 


4.37(0.25) 


3.6(0.3) 


67.76(0.15) 


3.15(0.28) 


- 


- 


- 


- 


?s 




MM3 


12.9(0.3) 


68.73(0.06) 


4.33(0.14) 


11.7(0.5) 


68.52(0.09) 


4.21(0.21) 


- 


- 


- 


- 


CTQ 




MM4 


8.9(0.4) 


68.90(0.09) 


4.42(0.20) 


7.2(0.5) 


68.84(0.13) 


4.26(0.32) 


- 


- 


- 


- 




G022.06+00.21.. 


MMl 


16.3(0.4) 


51.24(0.06) 


5.33(0.15) 


27.0(0.9) 


51.12(0.09) 


5.31(0.22) 


- 


- 


- 


- 


= • 


G023.28-00.12.. 


MMl 


4.7(0.4) 


98.96(0.14) 


4.01(0.40) 


5.1(0.6) 


98.57(0.18) 


3.03(0.46) 


2.01(0.51) 


97.8(0.2) 


- 


- 


G024.37-00.15.. 


MM2 


5.4(0.3) 


56.21(0.11) 


3.92(0.27) 


6.1(0.7) 


57.11(0.39) 


5.48(0.86) 


- 


- 


- 


- 




G024.61-00.33.. 


MMl 


6.7(0.3) 


42.86(0.07) 


3.49(0.19) 


9.8(0.6) 


42.80(0.10) 


3.10(0.24) 






3.08(0.54) 


43.5(0.2) 


3 
p. 


G024.94-00.15.. 


MMl 


7.4(0.3) 


A '-I 11 / r\ r\r\\ 

47.13(0.09) 


3.94(0.22) 


8.7(0.6) 


A '—ir\/r\ 1 ^\ 

46.70(0.17) 


4.0 / (,UOJ ) 




45.9(0.2) 


2.01(0.44) 


47.9(0.2) 


O 




MMl 


7 7C0 '\\ 


7S ^Of'O 07~> 




7 7^0 1\ 


7S 17f'0 \ X\ 


3.14(0.34) 










S- 


G030.90+00.00B. 


MM2 


6.7(0.3) 


92.96(0.11) 


4.69(0.26) 


5.2(0.6) 


92.53(0.19) 


3.54(0.54) 










o 

i-b 


G034.71-00.63.. 


MMl 


6.1(0.4) 


44.76(0.12) 


4.40(0.29) 


16.6(0.9) 


44.38(0.15) 


6.16(0.41) 


3.15(0.50) 


43.8(0.2) 


3.17(0.50) 


45.8(0.2) 






MM2 


2.7(0.4) 


45.88(0.29) 


4.09(0.62) 


1.2(0.5) 


45.14(0.13) 


0.74(0.44) 










3' 


G053. 8 1-00.00 


MMl 


3.7(0.34) 


24.44(0.11) 


2.48(0.32) 
















rmation. 



























Notes. Columns arc (from left to right) source name, N2H^(3-2) peak main beam brightness temperature, LSR velocity at the N2H^(3-2) peak intensity, FWHP HCO*(4-3) line width, N2H^(3-2g^ 
peak main beam brightness temperature, LSR velocity at the HCO^(4-3) peak intensity, FWHP HCO^(4-3) line width, blue HCO^(4-3) line intensity peak and its LRS velocity, and red HCO^(4-3S' 
line intensity peak and its LSR velocity. 2. 
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Table B.6. Star formation tracers (ii): K levels of CH3CN(5^). /T-level detections are marked as 'V, non-detections as 



Source name 







K levels 
1 2 


3 


r„,(CH3CN) 
(K) 


7'.o,(NH3) 
(K) 


N(CH3CN) 
(lO'^cm-2) 


no 1 4. (S^ on S7 


MMl 


V 


V 


V 


V 




1 8 1 

1 0. i 






MM? 


^/ 

V 


V 


- 






1 S 7 






MM? 


^/ 

V 


V 


V 






1 8 7 
1 o. / 






MM3 


,/ 

V 


V 


- 






20.0 






MM4 


V 


V 


- 






7n 1 






MM1 

iVliVl i 


^/ 

V 


V 








1 S ? 






MM3 


V 


V 


V 






1 5 7 






MM4 


V 


V 








16.6 




G022.06+00.21... 


MMl 


V 


V 


V 


V 


32(5) 


24.7 


9.7(2.8) 


G023.28-00.12... 


MMl 


V 


V 








18.1 




G024.37-00.15... 


MM2 


V 


V 








15.7 




G024.61-00.33... 


MMl 


V 


V 


V 


V 


68(28) 


17.5 


5.9(4.7) 


G024.94-00.15... 


MMl 


V 


V 


V 




39(31) 


15.2 






MM2 


V 


V 








15.2 




G034.71-00.63... 


MMl 


V 


V 


V 


V 


29(6) 


17.8 


5.4(2.8) 


G035.49-00.30A.. 


MMl 


V 


V 








13.2 





Notes. Columns are (from left to right) source name, detection of CH3CN K levels, CH3CN rotational temperature, NH3 rotational temperature, 
and CH3CN column density. 

Table B.7. Star formation tracers (i): SiO(2-l) and H2CO(4o3-3o4) based on Gaussian fits. Non detections are marked by see 
Table IB. ll for upper limits. 



Source name 




Tmb 
(K) 


SiO(2-l) 

vlsr 
(km s-') 


(km s-') 




Tmb 
(K) 


H2CO(4-3) 
Vlsr 
(km s-') 


(km s-') 


GOB, 


.28-00.34... 


MMl 


0.11(0.03) 


39.47(0.43) 


4.74(0.88) 
















GOB, 


.91-00.51... 


MMl 











.53(0. 


17) 


22, 


.94(0. 


,13) 


1.24(0.3) 


G014, 


.39-00.75 A.. 


MMl 








0, 


.80(0. 


30) 


17, 


.64(0. 


,14) 


1.01(0.28) 


G014, 


.63-00.57... 


MMl 


0.45(0.04) 


17.60(0.10) 


5.17(0.34) 


2, 


.41(0. 


10) 


17, 


.23(0. 


,05) 


3.62(0.11) 






MM2 


0.08(0.02) 


18.73(0.53) 


6.42(1.28) 





.85(0. 


15) 


18, 


.24(0. 


,12) 


2.71(0.37) 


G017, 


.19+00.81... 


MMl 


0.09(0.02) 


24.68(0.50) 


7.77(1.05 


0, 


.56(0. 


B) 


24, 


.96(0. 


,17) 


2.28(0.39) 






MM2 


0.79(0.11)° 


22.22(0.25Jl^ 


26.48(0.64) 


1 


.39(0. 


10) 


22, 


.46(0. 


,12) 


5.92(0.34) 






MM3 


0.76(0.11)° 


21.95(0.25)° 


19.42(0.60) 
















G018, 


.26-00.24... 


MMl 


0.36(0.02) 


66.88(0.16) 


9.18(0.47) 


0, 


.69(0. 


12) 


68, 


.02(0. 


,18) 


3.69(0.49) 






MM2 


0.09(0.02) 


68.01(0.64) 


11.14(2.35) 


0, 


.41(0. 


12) 


67, 


.74(0. 


,29) 


3.12(0.70) 






MM3 


0.21(0.02) 


68.50(0.27) 


10.06(0.85) 


0, 


.46(0. 


14) 


68, 


.52(0. 


,22) 


3.06(0.75) 






MM4 


0.25(0.02) 


69.24(0.28) 


13.02(0.83) 


0, 


.34(0. 


14) 


69 


.91(0. 


,43) 


4.30(1.46) 


G022, 


.06+00.21... 


MMl 


0.24(0.02) 


48.27(0.35) 


16.26(0.89) 


1 


.88(0. 


11) 


51 


.04(0. 


,08) 


4.78(0.22) 


G023, 


.28-00.12... 


MMl 


0.13(0.02) 


101.05(0.55) 


12.37(1.39) 
















G024, 


.37-00.15... 


MM2 


0.25(0.02) 


54.53(0.30) 


10.81(0.84) 


0, 


.69(0. 


11) 


55, 


.79(0. 


21) 


5.13(0.64) 


G024, 


.61-00.33... 


MMl 


0.11(0.02) 


40.44(0.49) 


10.18(1.2) 


0, 


.54(0. 


10) 


42, 


.96(0. 


24) 


3.63(0.50) 


G024, 


.94-00.15... 


MMl 


0.24(0.03) 


47.14(0.19) 


5.91(0.51) 


0, 


.69(0. 


12) 


46, 


.83(0. 


19) 


3.87(0.50) 






MM2 


0.18(0.02) 


47.39(0.38) 


11.61(1.08) 
















G025, 


.79+00.81... 


MMl 


0.12(0.04) 


49.20(0.56) 


8.04(2.01) 
















G034, 


.71-00.63... 


MMl 


0.15(0.02) 


44.40(0.46) 


14.14(1.19) 


1 


.49(0. 


13) 


44, 


.12(0. 


,08) 


3.02(0.21) 






MM2 


0.35(0.08) 


45.47(0.23) 


6.63(1.25) 


0, 


.27(0. 


09) 


45, 


.26(0. 


,57) 


5.03(1.31) 






MM3 


0.17(0.04) 


45.94(0.20) 


3.37(0.62) 
















G035, 


.49-00.30B.. 


MM2 


0.14(0.03) 


45.07(0.25) 


3.88(0.67) 




















MM3 


0.18(0.05) 


45.32(0.14) 


1.57(0.32) 
















G053, 


.81-00.00... 


MMl 


0.11(0.01) 


23.14(0.24) 


9.14(0.52) 


0, 


.47(0. 


11) 


23, 


.79(0. 


,39) 


4.81(0.79) 



Notes. Columns are (from left to right) source name, SiO peak main beam brightness temperature, LSR velocity at the SiO peak intensity, FWHP 
SiO line width, H2CO peak main beam brightness temperature, LSR velocity at the H2CO peak intensity, and FWHP H2CO line width. 
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Table B.8. Observed properties of J8o(2-l) based on Gaussian fits. Non detections are marked by see Table IB. ll for upper limits. 



Source name 








Av 








(Kkms"') 


(km s-') 


(kms-i) 




7Q nn '^s 


MMl 


- 


- 


- 






MM2 


2.1(0.7) 


6.14(0.44) 


1.84(0.7) 




nn 'XA 


MMl 


4.3(0.5) 


41.09(0.12) 


1.53(1.79) 




Q 1 nn 1 


MMl 


9.7(0.5) 


22.95(0.07) 


2.81(0.19) 




07 nn 


MMl 


12.4(0.5) 


19.66(0.05) 


2.57(0.12) 


nni A 


■iQ nn 7^ A 

.jy—y}\j. /JJ\.. 


MMl 


12.9(0.5) 


17.67(0.04) 


2.10(0.10) 






MM2 


14.1(0.5) 


17.6(0.04) 


2.15(0.08) 


nni A 


nn 7^Ti 


MM3 


5.8(0.3) 


21.06(0.04) 


1.81(0.28) 


nni A 


nn '^7 


MMl 


18.8(0.2) 


18.02(0.02) 


3.14(0.05) 






MM2 


17.1(0.3) 


18.34(0.02) 


2.77(0.05) 






MM3 


11.0(0.2) 


17.36(0.02) 


2.24(0.04) 






MM4 


6.2(0.2) 


18.9(0.03) 


1.96(0.14) 


nni fx 


.yJ•T\J^J.z.^■... 


MMl 


8.7(0.3) 


23.55(0.04) 


1.99(0.08) 


nni 7 


.ly+uu.oi... 


MMl 


9.9(0.2) 


24.23(0.02) 


2.85(0.06) 






MM2 


9.6(0.2) 


22.63(0.02) 


2.41(0.05) 






MM3 


9.0(0.1) 


22.75(0.02) 


2.51(0.05) 






MM4 


12.5(0.1) 


21.72(0.01) 


2.66(0.04) 


nni s 


OA nn '7,1 


MMl 


10.3(0.5) 


67.82(0.09) 


3.70(0.25) 






MM2 


18.2(0.4) 


67.44(0.03) 


2.86(0.08) 






MM3 


15.5(0.4) 


68.51(0.05) 


3.72(0.11) 






MM4 


13.3(0.7) 


68.4(0.1) 


3.78(0.22) 






MM5 


17.0(0.5) 


66.48(0.04) 


2.82(0.08) 


nmo 


.UO+UU.Zl... 


MMl 


19.0(0.3) 


51.29(0.02) 


3.08(0.05) 






MM2 


9.6(0.4) 


51.87(0.05) 


2.69(0.12) 




o Q nn 1 1 
.Zo— uu.iz... 


MMl 


14.7(0.9) 


99.15(0.08) 


2.89(0.22) 






MM2 


6.1(0.8) 


99.54(0.15) 


1.96(0.27) 




T7 nn 1 ^ 


MMl 


14.6(0.4) 


58.99(0.04) 


2.93(0.1) 






MM2 


4.6(0.5) 


56.32(0.2) 


3.67(0.44) 




Ai nn 
.01— UU.jj... 


MMl 


12.9(0.2) 


42.43(0.02) 


2.45(0.04) 






MM2 


3.9(0.3) 


43.43(0.09) 


3.12(0.26) 


nnoA 


QA nn 1 ^ 


MMl 


5.5(0.3) 


47.07(0.05) 


2.48(0.16) 






MM2 


4.7(0.3) 


47.74(0.1) 


3.26(0.2) 


nno'N 
LrUZj, 


7Q_i_nn G 1 


MMl 


3.4(0.2) 


49.8(0.06) 


2.16(0.19) 






MM2 


4.2(0.2) 


49.04(0.05) 


2.77(0.12) 


nmn 


Qn_i_nn nn a 


MMl 


19.8(0.7) 


75.11(0.05) 


3.22(0.13) 


nn^n 


Qnj_nn nnR 


MM2 


17.8(0.4) 


92.64(0.04) 


3.86(0.10) 


nn^n 


Qn_i_nn nnn 


MM3 


15.4(0.5) 


94.35(0.08) 


4.46(0.18) 




71 nn A^ 

. / 1— UU.Oj... 


MMl 


7.9(0.2) 


44.93(0.04) 


3.31(0.12) 






MM2 


5.7(0.2) 


45.96(0.06) 


2.82(0.14) 






MM3 


6.3(0.2) 


46.55(0.04) 


2.66(0.09) 


G034, 


.77-00.81... 


MMl 


5.2(0.3) 


43.72(0.08) 


3.44(0.21) 


G034, 


.85+00.43... 


MMl 


2.8(0.2) 


55.51(0.07) 


2.11(0.14) 


G035, 


.49-00.30A.. 


iVliVl 1 




zLf'n n^^ 


1 lliC\ C\f\\ 


G035, 


.49-00.30B.. 


MM2 


4.1(0.2) 


45.16(0.04) 


2.04(0.13) 






MM3 


5.7(0.2) 


45.67(0.05) 


2.46(0.11) 


G037, 


.44+00. 14A.. 


MMl 


1.0(0.1) 


18.07(0.05) 


1.37(1.5) 


G037, 


.44+00. 14B.. 


MM2 


8.0(0.2) 


40.11(0.03) 


1.87(0.05) 


G050 


.06+00.06... 


MMl 


10.1(0.2) 


54.19(0.02) 


2.27(0.05) 






MM2 


10.4(0.3) 


54.82(0.03) 


2.06(0.05) 


G053, 


.81-00.00... 


MMl 


3.8(0.1) 


23.89(0.02) 


2.30(0.07) 



Notes. Column are (from left to right) source name, integrated line intensity, LSR velocity at the line intensity peak, and FWHP line width. 
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Table B.9. Derived abundance and depletion from the 18o(2-l) line. 



Source name 




T'mlCNH,) 






Vr 1 8 (\ 


rj 






(K) 


( lO'" cur-) 


(10" cnr-) 


( 1 0"^ ) 




G012.73-00.58... 


MMl 


9.3 


- 


4.5 


- 


- 




MM2 


11.4 


1.3(0.5) 


- 


- 


- 


G013.28-00.34... 


MMl 


14.9 


2.4(0.3) 


2.9 


8.3 


2.1 


G013.91-00.51... 


MMl 


14.1 


5.5(0.3) 


1.1 


7.1 


2.4 


G013.97-00.45... 


MMl 


16.6 


6.9(0.3) 


- 


- 


- 


G014.39-00.75A.. 


MMl 


17.4 


7.2(0.3) 


5.7 


12.6 


1.3 




MM2 


20.0 


7.9(0.3) 


- 


- 


- 


G014.39-00.75B.. 


MM3 


11.8 


3.5(0.2) 


5.7 


6.2 


2.8 


G014.63-00.57... 


MMl 


18.1 


10.5(0.1) 


26.3 


4 


4.3 




MM2 


15.7 


9.6(0.2) 


15.4 


6.2 


2.7 




MM3 


15.8 


6.1(0.1) 


5.2 


11.8 


1.4 




MM4 


19.1 


3.4(0.1) 


3 


11.5 


1.5 


G016.93+00.24... 


MMl 


14.0 


4.9(0.2) 


4.2 


11.8 


1.4 


G017.19+00.81... 


MMl 


17.2 


5.5(0.1) 


4.9 


11.3 


1.5 




MM2 


18.7 


5.4(0.1) 


17.4 


3.1 


5.5 




MM3 


20.0 


5.1(0.1) 


5.3 


9.6 


1.8 




MM4 


20.1 


7.0(0.1) 


2.9 


24.3 


0.7 


G018.26-00.24... 


MMl 


18.2 


5.7(0.3) 


9.8 


5.9 


2.9 




MM2 


17.4 


10.1(0.2) 


8.9 


11.4 


1.5 




MM3 


15.7 


8.7(0.2) 


7.1 


12.2 


1.4 




MM4 


16.6 


7.4(0.4) 


5.7 


12.9 


1.3 




MM5 


16.8 


9.5(0.2) 


6 


15.8 


1.1 


G022.06+00.21... 


MMl 


24.7 


11.2(0.2) 


25.4 


4.4 


3.8 




MM2 


15.5 


5.4(0.2) 


7.9 


6.8 


2.5 


G023.28-00.12... 


MMl 


18.1 


8.2(0.5) 


6.5 


12.6 


1.3 




MM2 


10.0 


3.4(0.6) 


- 


- 


- 


G024.37-00.15... 


MMl 


18.6 


8.1(0.2) 


4.9 


16.6 


1 




MM2 


15.7 


2.6(0.3) 


4.4 


5.8 


2.9 


G024.61-00.33... 


MMl 


17.5 


7.2(0.1) 


5.9 


12.2 


1.4 




MM2 


15.5 


2.2(0.1) 


3.5 


6.3 


2.7 


G024.94-00.15... 


MMl 


15.2 


3.1(0.1) 


6.5 


4.7 


3.6 




MM2 


15.2 


2.7(0.1) 


5.5 


4.8 


3.5 


G025.79+00.81... 


MMl 


13.1 


1.9(0.1) 


- 


- 


- 




MM2 


15.0 


2.4(0.1) 


- 


- 


- 


G030.90+00.00A.. 


MMl 


18.6 


11(0.4) 


5.2 


21.2 


0.8 


G030.90+00.00B.. 


MM2 


15.0 


10(0.2) 


- 


- 


- 


G030.90+00.00C.. 


MM3 


16.4 


8.6(0.3) 


4.3 


20 


0.9 


G034.71-00.63... 


MMl 


17.8 


4.4(0.1) 


5.8 


7.6 


2.2 




MM2 


12.4 


3.4(0.1) 


5.2 


6.5 


2.6 




MM3 


17.1 


3.5(0.1) 


2.5 


14 


1.2 


G034.77-00.81... 


MMl 


15.0 


2.9(0.1) 


- 


- 


- 


G034.85+00.43... 


MMl 


13.2 


1.6(0.1) 










IVllVll 


io.O 


f, f,(n 1 

O.Ol^U. 1 ) 


o.o 


in 


1 7 
1. / 


G035.49-00.30B.. 


MM2 


11.9 


2.5(0.1) 


5.3 


4.6 


3.7 




MM3 


13.6 


3.2(0.1) 


4.2 


1.1 


2.2 


G037.44+00.14A.. 


MMl 


15.0 


0.6(0.05) 








G037.44+00.14B.. 


MM2 


15.0 


4.5(0.1) 








G050.06+00.06... 


MMl 


14.9 


5.6(0.1) 


4.4 


12.8 


1.3 




MM2 


14.1 


5.9(0.1) 


4.6 


12.9 


1.3 


G053.81-00.00... 


MMl 


12.4 


2.2(0.1) 


6.9 


3.2 


5.3 



Notes. Columns are (from left to right) source name, NH3 rotational temperature, J 80 column density, hydrogen column density from the 1.2 mm 
emission (Paper I), }8o abundance, and CO depletion factor. 



